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I.  INTRODUCTION 


A  guided  missile  is  one  which  is  usually  fired  in  a  direction 
approximately  towards  the  target  and  subsequently  receives  steering 
commands  from  the  guidance  system  to  improve  its  accuracy. 

In  recent  years,  the  application  of  Bank-to-Tum  (BTT)  guidance  to 
tactical  missiles  has  generated  considerable  interest,  motivated  by  certain 
advantages  that  a  Bank-to-Tum  control  configuration  can  offer. 

First,  against  high  performance  threats,  there  is  a  need  for  defensive 
missiles  to  develops  increasingly  higher  lift  accelerations,  which  requires 
banking  maneuver  to  properly  direct  the  control  vector. 

Second,  the  need  for  greatly  increased  stand  off  ranges  has  led  the 
design  of  tactical  missiles  towards  airbreathing  propulsion  systems  such  as 
the  ramjet.  These  designs  generally  have  configuration  geometries  that  are 
not  cruciform,  due  to  exposed  inlets  beneath  the  vehicle.  As  a  result,  there 
are  often  stringent  limits  on  the  sideslip  angle  that  can  be  developed  during 
engine  operation,  and  this  too,  dictates  some  types  of  Bank-to-Tum  control 
scheme. 

This  thesis  investigates  the  controllability  of  a  Bank-to-Tum  missile, 
when  one  or  more  control  surfaces  are  damaged,  at  fixed  flight  conditions 
(i.e.  constant  altitude,  Mach  number  and  missile  weight  and  inertias). 

The  approach  is  to  design  a  control  reconfiguration  module  to  be  placed 
between  the  flight  control  laws  and  the  control  surfaces.  In  the  face  of 
missile  damage  the  reconfiguration  module  would  reallocate  control 


commands  to  the  unimpaired  control  surfaces  so  that  flight  critical  pitch, 
yaw,  and  rolling  moments  would  he  preserved  to  the  maximum  extent 
possible.  Using  the  reconfiguration  module  control  laws,  the  existing  flight 
control  laws  would  not  be  altered. 

In  the  first  part  of  the  thesis  the  control  surfaces  conventions, 
definitions  and  the  manoeuvre  by  banking  are  reviewed. 

The  second  and  third  part  reviews  the  design  of  the  longitudinal  and 
lateral  channels  and  involves  the  uncoupled  pitch,  roll  and  yaw  channels  for 
the  circular  and  elliptical  airframes  of  a  B.T.T.  missile.  The  resultant 
continuous  open  loop  designs  are  analyzed  and  found  to  be  in  accordance 
with  the  desired  specified  requirements.  Applying  the  convention  control 
theory,  the  corresponding  modified  longitudinal  and  lateral  autopilots  ore 
obtained,  analyzed  and  compared  with  the  original  system.  The  effect  of 
damage  to  one  or  more  control  surfaces,  on  the  controllability  of  each 
system,  is  obtained,  analyzed  and  compared  with  the  unimpaired  models 
respectively. 
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Finally,  the  combined  linear  open  loop  pitch  and  roll  channel  autopilots, 
for  the  circular  airframe  configuration,  ore  Investigated.  The  effect  of 
damage  to  one  or  more  actuators,  on  the  control  of  the  combined  unimpaired 
system  is  obtained  and  analyzed.  The  most  severe  type  of  impairment 
considered  In  this  thesis,  l.e.,  the  hard  over  type.  Is  extensively  investigated 
in  Chapter  V.  Corrective  actions  are  then  developed  and  presented. 

The  robustness  of  the  complete  system  ts  analyzed  in  chapter  five. 
Conclusions  and  recommendations  for  future  study  are  then  given. 
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The  Investigation  and  analysis  In  all  the  above  cases  was  performed 
using  the  existing  at  Naval  Postgraduate  school  (NPS),  optimal  Systems 
Control  Program  (OPTSYS)  for  the  continuous  system.  Robustness  Design 
Program  (POPLAR)  for  the  robustness  of  the  system,  and  EASYPLOT 
program  for  the  plotting  of  different  kind  of  Impairments  and  comparisons. 
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DU  D  A  ill/  |||^ 

BY  BANKING 

In  this  chapter,  before  going  Into  mathematical  details  concerning  the 
motion  of  a  guided  missile  in  space  as  a  result  of  guidance  commands,  some 
definitions  and  discussion  concerning  the  control  surfaces  are  given. 

A.  GENERAL 

The  task  of  a  Missile  Control  system  Is  to  manoeuvere  the  missile 
quickly  and  efficiently  as  a  result  of  the  following  two  control  slngnals: 

(a)  Detection  whether  the  missile  is  flying  too  high  or  too  low,  or  too 
much  to  the  left  or  right. 

(b)  Measure  the  deviations  or  errors  and  develope  commands  to  the 
control  system  to  reduce  these  errors  to  zero. 

B.  CONTROL  SURFACES  DEFINITIONS 

Suppose  that  the  control  surfaces  1  and  3,  as  Is  shown  in  the  Figure 
2.1,  are  mechanichally  linked  together  such  that  a  servo  must  impart  the 
same  rotation  to  both.  Then  these  surfaces  are  essentially  elevators. 
Positive  elevator  deflection  produces  a  negative  force  in  the 
z-dl  recti  on  and  an  anti-clockwise  moment  about  the  y-axls. 

If  now  surfaces  1  and  3  each  have  their  own  servo  it  is  possible  for 
them  to  act  as  ailerons.  Positive  aileron  deflection  produces  an  anti¬ 
clockwise  moment  about  the  x-axis. 
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If  looking  in  the  direction  y  one  surface  Is  rotated  6°  clockwise  and 
the  other  surface  6°  anti-clockwise  then  a  pure  couple  is  Imparted  to  the 
missile  about  the  fore  and  aft  axis  and  this  will  tend  to  moke  the  missile 
roll.  Such  control  surfaces  are  now  called  ailerons,  and  by  doing  the  same 
thing  to  control  surfaces  2  and  4,  the  turning  effect  of  the  ailerons  con  be 


doubled. 
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Figure  2.1  Control  surfaces  looking  from  rear  of  the  Missile. 

If  the  aerodynamics  are  linear  i.e.,  the  normal  forces  are  proportional 
to  incidence  then  the  principle  of  superposition  applies.  Commmands  for 
elevator  rudder  and  aileron  movements  can  be  added  electrically  resulting  In 
unequal  movements  to  opposite  control  surfaces.  Positive  rudder 
deflection  produces  a  positive  force  In  the  y-dlrection  and  negative 
moment  about  the  z-axis.  In  this  way  we  have  the  means  to  control  roll 
motion  as  well  as  the  up-down  i.e.,  pitch  motion  and  left-right  i.e.,  yaw 
motion. 


C.  CONTROL  SURFACE  CONVETION 

As  is  shown  in  the  Figure  2.1  deflections  6t,  62,  63,  64  are  defined  as 
positive  If  the  leading  edge  of  the  control  surface  is  rotated  clockwise  as 
one  looks  outwards  along  the  Individual  hinge  axis.  The  total  aileron 
deflection  is  defined  as: 

Aileron  deflection  6=  1 /4  (6,  +62+63+54)  (II.C- 1 ) 

or  1/2  (6, +63)  or  1/2(62+64)  (II.C-2) 

if  only  two  surfaces  act  differentially. 

D.  BANK-TO-TURN  CONTROL  SURFACES  CONFIGURATION 

The  same  kind  of  analysis  can  be  applied  for  the  Bank-to-Tum  missile 
configuration,  shown  In  Figure  2.2.  Note  that  the  arrow  on  the  control 
surface  indicates  the  direction  of  positive  leading  edge  deflection. 


Figure  2.2  B.T.T  Aerodynamic  Sign  Convention. 


The  maneuver  plane  Is  the  x8-zg  plane.  The  downwards  movement  is 
defined  as  positive.  This  is  the  configuration  investigated  in  this  thesis. 


The  B.T.T.  missile  configuration  Is  tail-controlled  using  four  Identical 
control  surfaces.  The  control  surfaces  are  located  flush  with  the  body  base 
with  a  ±30  degrees  dihedral.  With  tall  control  configurations  the  pitching 
moment  or  control  effectiveness  Is  quite  high  In  view  of  the  fact  that  the 
tall  moment  arm  is  generally  fairly  large.  Also  there  is  no  adverse 
downwash  from  deflection  of  the  forward  surfaces.  In  comparison  with  the 
canard  control  configuration.  A  more  complete  model  of  the  B.T.T.  circular 
configuration  Is  also  shown  In  Figures  3.1  and  3.4 ,  taken  from  Arrow,  A. 
iRef.  11. 

E.  MANEUVER  BY  BANKING. 

The  Conventional  method  of  altering  course  to  the  right  or  left  in  a 
missile  or  aircraft  is  to  use  the  ailerons  to  bonk  l.e,  roll  by  an  angle  <p,  as  Is 
shown  in  Figure  2.3.  If  the  lift  force  Is  Increased  slightly  by  the  use  of  the 
elevators  so  that  the  vertical  component  of  lift  equals  the  weight  then 
there  Is  a  horizontal  component  of  lift  equal  to  the  total  lift  times  sin  6.  It 
Is  this  component  of  lift  which  causes  the  flight  path  to  change.  An  exact 
analysis  of  the  maneuver  Is  not  simple  as  the  airflow  over  the  outer  control 
surface  is  faster  than  that  over  the  inner  control  surface  due  to  the  rate  of 
turn.  In  practice  a  small  amount  of  rudder  Is  applied  In  an  attempt  to  make 
the  general  airflow  directly  along  the  fore  and  aft  axis  of  the  cordinated 
B.T.T.  missile  and  In  the  plane  of  the  control  surfaces.  In  this  condition 
there  Is  no  ’sideslip’  and  hence  no  net  sldeforce.  This  Is  the  preferred 
method  of  manoeuvring  since  lifting  forces  are  most  efficiently  generated 
perpendicular  to  the  control  surfaces;  the  llft-to-drog  ratio  Is  a  maximum 


In  this  condition.  Even  so  there  Is  some  additional  Induced  drag  when 
manoeuvring.  An  additional  manoeuvring  force  is  obtained  by  increasing  the 
Incidence  and  this  means  that  the  normal  force  Is  Inclined  backwards  to  the 
velocity  vector  by  the  angle  of  incidence.  There  Is  a  useful  component  equal 
to  the  normal  force  times  the  cosine  of  the  angle  of  incidence  and  a 
component  proportional  to  the  sine  of  the  incidence  opposing  the  velocity 
vector.  This  ‘Induced  drag*  will  always  occur  whenever  a  lifting  force  Is 
produced  irrespective  of  the  method  of  manoeuvre. 


Figure  2.3  Forces  Acting  on  a  Banking  Missile  and  Coordinated  Missile 
Motion  for  coordinated  B.T.T.  Control  Policies. 

Becouse  it  Is  desired  that  the  B.T.T.  missile  avoid  negative  angles  of 
attack,  the  missile  is  forced  to  roll  about  its  velocity  vector  when  the 
desired  manouvre  direction  Is  In  the  negative  angle-of -attack  direction. 
Although  some  negative  angles  of  attack  may  be  possible  with  certain 
missiles. 


Coordinated  motion  or  zero  sideslip  angle  Is  achieved  by  directing  the 
body  fixed  pitch  axis  of  the  missile  at  the  missile  velocity  vector  so  that 
there  Is  no  component  of  missile  velocity  vector  along  the  body  fixed  yaw 
axis  of  the  missile.  The  Figure  2.3  (b)  shows  the  attitude  of  the  missile 
body  with  respect  to  Its  velocity  vector  v?  When  an  upward  command  is 
applied,  f.e.,  <p  =  0,  the  missile  body  moves  upward  with  Its  pitch  axis 
directed  at  the  velocity  vector  until  It  reaches  manoeuver  level  or  angle  of 
attack.  No  roll  motion  Is  required  to  maintain  coordination  of  this 
manoeuver.  The  circular  airframe  configuration  of  this  thesis  Is  a  planar 
configuration  which  Is  symmetric  about  Its  wing  planes  and  therefore  have 
two  perfect  directions  which  are  normal  to  the  planes  of  its  wing. 
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A.  INTRODUCTION 

In  this  Chapter  the  design  of  the  continuous  uncoupled  pitch  channel 
(longitudinal)  autopilot  Is  reviewed,  as  was  analyzed  by  Arrow,  A.  and 
Karalskos,  I.  [Ref.  1  and  21.  A  modified  longitudinal  autopilot  Is  then 
designed  with  the  four  control  surfaces  split  Into  independent  components, 
for  the  circular  configuration  model  shown  In  Figure  3.1. 

Then  a  comparison  between  the  two  systems  (original  and  modified) 
shows  that  are  Identical.  The  transient  response  plots  and  their  analysis 
demonstrate  the  performance  of  the  modified  pitch  channel. 

Having  established  the  validity  of  the  modified  system  one  then  uses 
the  multiple-input  multiple-output  control  theory  and  state-space 
representation  [Ref.  3],  to  investigate  the  performance  of  the  modified  pitch 
channel  autopilot  for  various  classes  of  damages. 

For  this  thesis,  two  classes  of  Impairments  are  considered  namely, 
locked  or  hard  over  and  inoperative  with  the  control  surfaces  centered.  The 
plots  of  a  modified  pitch  autopilot  without  Impairment  are  used  as  a  basis 
for  comparison,  the  Impairment  induced  effects  on  normal  acceleration,  rate 
and  tail  Incidence  are  assessed. 

A  general  block  diagram,  of  a  B.T.T.  autopilot  Is  shown  in  Figure  3.2 
taken  from  [Ref.  1].  Inertial  acceleration  commands  are  applied  In  polar 
coordinates  (l.e.,  magnitude  of  the  command  (  %  )  applied  to  the  pitch 
autopilot  and  the  direction  (  fc )  Is  applied  to  the  roll  autopilot.  Achieved 


manoeuver  plane  or  Inertial  acceleration  In  rectangular  coordinates  ( i.e.,  t\z 
and  i\v)  is  determined  by  resolving  achieved  body-fixed  accelerations  ( f.e., 
Hzond  T)v)  through  missile  roll  angle  (f). 

Before  going  into  the  details  concerning  the  various  kinds  of  damages 
It  Is  prefered  to  make  here  clear  some  terms  used  in  this  work. 

ORIGINAL  SYSTEM  is  called  the  unmodified  pitch  channel 
autopilot  presented  in  (Ref.  1]. 

UNIMPAIRED  SYSTEM  is  called  the  modifiedy  i.e.,  with  the  split 
control  surfaces  system,  and  without  any  kind  of  damage. 

IMPAIRMENT,  HARD  OVER  OR  INOPERATIVE  are  terms  expressing 
a  specific  kind  of  damage  to  one  or  more  control  surfaces. 

B.  UNCOUPLED  LINEAR  PITCH  CHANNEL  (ORIGINAL) 


The  linear  pitch  channel  dynamic  model  which  is  used  for  the 
circular  airframe.  Is  shown  In  Figure  3.3,  with  the  following  requirements : 

a.  In  actuator  command  branch,  the  high  attenuation  frequency  must  be  i 
15  db  at  100  rad/sec  and  zero  angle-of-ottack.  T.tis  requirement 
limits  the  autopilot  speeds  of  time  response. 

b.  Relative  stability  with  gain  margins  i  6  db,  phase  margins  i  30 
degrees  with  a  goal  at  12  db  and  50  degrees. 

c.  An  acceleration  time  responce  with  63X  time  constant  of  0.5  seconds 
for  a  step  command  of  acceleration  at  the  flight  condition  of  interest, 
overshoot  i  10%,  and  zero  steady-state  error  to  reduce  variations  of 
guidance  navigation  gain. 
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The  linear  pitch  channel  dynamic  model  which  is  used,  is  shown  in 
Figure  3.3.  This  Figure  is  token  from  [Ref.  1].  The  model  is  linearized  for 
stability  studies  and  the  following  three  assumptions  were  made : 

a.  Manoeuver  plane  Is  the  plane  xg  -  Ze  ,  as  Is  shown  in  Figure  3.4  in 
which  the  arrow  indicates  the  direction  of  positive  leading  edge 
deflection. 

b.  The  missile  is  trimmed  in  pitch,  at  fixed  values  of  angle  of  attack  (a), 
pitch  angular  rate  (q),  and  pitch  tail  incidence  (6p),  that  is  0. 

c.  The  roll  rote  (P)  of  the  missile  is  constant. 

Using  the  foregoing  assumptions,  and  the  approximation  that  the 
variation  In  forward  velocity  (u),  is  equal  to  zero  a  simplified  set  of 
longitudlnol  equations  of  motion  ore  constructed  oppllcoble  to  the  short 
period  of  motion.  These  equations  describe  the  two  degree  of  freedom,  short 
period  longitudinal  mode. 

Linearized  aerodynamic  derivatives  for  M  =  3.95  end  Altitude  =  60 
KFT  are  provided  in  Toble  3.1,  as  is  given  in  [Ref.  1  ]. 

2.  pitch  Control  Low 

The  longitudinal  control  law  model  for  the  circular  airframe  was 
developed  from  the  block  diagrams  and  information  provided  in  [Ref.  1],  and 
is  shown  in  Figure  3.5.  Lead-log  filters  were  used  in  the  system  to  shape 
the  output  response,  provide  adequate  gain  and  phose  margins,  ond  finally  to 
prevent  guidance  noise  saturation  problems.  The  guidance  requirement  of 
zero  steady-state  error  is  satisfied  by  using  an  integrator  in  the 
acceleration  error  branch  of  the  circulore  control  law.  The  rate  error 
conpensatlon  determines  the  high  frequency  attenuation  ond  was  used  to 
minimize  the  effect  of  aerodynamic  variations  on  acceleration  time 


vs 


response.  The  acceleration  error  compensation  determines  the  acceleration 
time  response. 


The  blocks  In  the  pitch  control  law  of  Figure  3.5  for  the  circular 
model  represent  the  basic  transfer  functions  In  which  a  normal  acceleration 
command  (n&j)  equal  to  1  gee,  Is  applied  to  the  pitch  control  low  which  uses 

measurements  of: 

a.  Missile  body  pitch  angular  rate  (q)and 

b.  Pitch  normal  acceleration  (rfc) 

to  determine  the  required  actuator  comend  (6pc).  The  actuator  is  modeled  as 

a  first  order  lag  at  188.4  rad/sec  and  Is  shown  in  Figure  3.2. 

3.  Basic  Transfer  functions  of  the  Aerodynamic  Model 

The  equations  which  are  represented  by  the  block  diagram  of  the 
aerodynamic  model,  of  the  uncoupled  pitch  channel  are  presented  below  and 
taken  from  lRef.2l . 

a.  The  pitch  aerodynamic  transfer  function  for  the  pitch 
ongular  rote  (q)  about  the  Ye  axis: 

q  -1.3361s-  0.15 

—  - (1/sec)  (1 1 1.8.3- 1) 

6p  22.545  X  10"3  s*  ♦  3.3634  X  10*  3  ♦  1 

b.  The  pitch  aerodynamic  transfer  function  for  the  achieved 
moneuver  normal  acceleration  (r^)  in  the  Zq  direction: 

n,  1 06.074  x  10"*  s2-  18.82 

- - - (gs/rad)  (III.B.3-2) 

5p  22.545  X  10“5  s2  ♦  3.3634  X  10"5  s  ♦  1 


Figure  3. 1  B.T.T.  Model  of  Circular  Configuration;  Scale  1/6. 
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The  analysis  of  the  uncoupled  pitch  channel  autopilot  is  based  on 
the  time  responses  and  Bode  plots  of  maneuver  plane  acceleration,  body 
angular  rates  and  tall  Incidence  angles,  obtained  by  OPTSVS  control  program 
of  Naval  Postgraduate  School  (NPS),  which  requires  the  differential 
equations  to  be  In  the  state-space  variable  format: 


X  =  F  X  ♦  G  U 


(III.B.4- 1 )  to 


define  the  system.  In  this  system  the  F  matrix  is  the  open  loop  dynamics 
matrix  (system  or  plant)  and  the  6  matrix  is  the  control  matrix.  The 
variable  assignements  are  X  as  the  state  vector  and  u  as  the  control  input 
vector.  It  follows  that  X  Is  the  time  derivative  of  X. 

Using  state-space  representation  of  a  system,  in  which  the 
forcing  function  Involves  derivatives,  the  following  equations  are  derived 
and  can  be  modeled  as  a  tenth  order  system,  for  the  uncoupled  pitch 
channel : 


X,=  -150  X,  ♦ISO  Xfi  ♦  705.75  5p 

(III.B.4-4) 

X2  =  X3  -  1.3531  nzo  ♦  1.3531  X, 

(III.B.4-5) 

X3  =  -5  X3  ♦  6.572  nzc  -  6.572  X, 

(III.B.4-6) 

X,  =  Xg  -  59.2637  8p 

(1  MB  4-7) 

Xg  =  -M9. 186  x  10-3  X5  -44.3557  X4  +  1.7887  5p 

(III.B.4-0) 

Xg  =  X7  -  701.92*10-3  6p 

(III  B  4-9) 

X7  =  -149.106x10-3  X7  -  44.3557  Xg  -  1043.36  8p 

(III.B.4- 10) 

'$p  =  -188.4  5p  ♦  188.4  6po 

(III.B.4- 11) 

$Po  =  -0.M3  fipo  -2.2308  X2  ♦  2.2308  X4  -  340.58x10-3  X3 
♦  460.839x10-3  na,  -  460.839x10-3  X,  ♦  340.58 
x  1 0-3  Xa  -  20. 1 84  5p  (III.B.4- 1 2) 

=  -(  X7  -  887. 1 24  6p  ♦  886.422  Sp,  )  (III.B.4- 1 3) 


where  the  state  variables  are: 

•  Xt  j  X2 .  pitch  angular  rate  state  variables 

•  *6,  h :  achlved  body  acceleration  state  variables 

•  X5  .  output  of  acceleration  filter 

•  Spo  :  input  command  In  actuator  network 

•  Sp  .  tall  Incidence  angle 

Expressing  the  above  equations  as  a  set  of  first  order  matrix 
in  the  format  of  (III.B.4- 1)  the  plant  and  control  matrices  are: 


-150  00  0  0  0  0  0  0  150 

1.353  01  0  0  0  0  0  00 

-6.572  0  -5  000  0  0  00 

0  0  0  0  1  0  0  -59.26  0  0 

F=  0  0  0  -44.36  -0.15  0  0  1.789  0  0 

0  0  0  0  0  0  1  -0.7019  0  0 

0  0  0  0  0  -44.36  -0.1492  -1043  0  0 

0  0  0  0  0  0  0  -188.4  188.4  0 

-0.461-2.231-0.3406  2.231  0.3406  0  0  -20.18  -0.143  0 

0  0  0  0  0  0  -1  887.1  -886.4  0 

k  « 

(III.B.4- 14) 


The  control  matrix  G  is  following: 


-1.3531 
6.572 
0 

6s  |  0  |  (III.B.4-15) 

0 
0 
0 

0.4608 


where:  U  =  rfe,  as  control  input  to  the  system.  (III.B.4-16) 

According  to  the  [Ref.  1],  using  a  step  Input  function  which 
represents  the  '1  gee  Command*  and  a  trim  angle  of  attack  a»=0,  the  time 

response  plots  of  the  state-variables  of  Interest  were  obtained  from  the 
OPTSYS  program  and  are  represented  at  the  Figures  3.6,  3.7  and  3.8.  Typical 
performance  criteria  used  to  characterize  the  transient  response  to  the  one 
gee  unit  step  input  Include  overshoot,  delay  time,  rise  time,  and  settling 
time.  The  above-mentioned  criteria  are  defined  with  respect  to  the  step 
response: 

•  ttaxfmum  Overshot.  The  maximum  overshoot  is  defined  os  the  largest 
deviation  of  the  output  over  the  step  Input  during  the  transient  state. 
The  amount  of  maximum  overshoot  Is  also  used  as  a  measure  of  the 
relative  stability  of  the  system.  The  moximum  overshoot  is 
represented  in  this  thesis  as  a  percentage  of  the  final  value  of  the 
step  response:  that  is, 

max.  overshoot 

percent  maximum  overshoot  - x  1 00%  (III.B.4- 1 7) 

final  value 


•  Delay  time.  The  delay  lime  Is  defined  as  the  time  required  for  the  step 
response  to  reach  50%  of  its  final  value. 

•  Rise  time  The  rise  time  is  defined  as  the  time  required  for  the  step 
response  to  rise  from  10%  to  90%  of  Its  final  value.  Frequently  the 
rise  time  Is  also  refered  to  as  the  time  constant  of  the  system. 

•  Settling  time  The  settling  time  Is  defined  as  the  time  required  for  the 
step  response  to  decrease  and  stay  within  a  specified  percentage  of 
Its  final  value.  For  this  thesis  Is  specified  arbitrarily  at  0.0 1%. 

The  four  quantities  defined  above  give  a  direct  measure  of  the 

transient  characteristics  of  the  system  to  the  step  Input: 

a.  The  pitch  normal  acceleration  (rfe  no  Nz),  Figure  3.6,  g's  vs. 

time  which  has: 

(1)  A  0.5  rise  time  or  time  constant. 

(2)  A  minor  0.005  steady  state  error  and  2.2%  overshoot. 

(3)  A  0.28  sec  delay  time  and  2  sec  settling  time. 

b.  The  pitch  angular  rate  (q  or  Q).  Figure  3.7.  rad/sec  vs.  time. 

which  has: 

(1)  The  maximum  body  angular  rate  (-0.17  rad/sec)  occurs  at  0.28  second 
when  the  tall  Incidence  Is  at  0.02  rad  angle  of  deflection. 

(2)  From  the  moment  at  0.28  sec  the  rate  Increases  and  reaches  at  2  sec 
Its  steady-state  value  (-0.01  rad/sec). 

(3)  The  -0.01  rad/sec  steady-state  constant  rate  Is  due  to  the  fact  that 
one  has  a  type  zero  system  ami  therefore  has  a  finite  error  to  o  step 
Input. 

c.  The  pitch  tail  incidence  (6p  or  DP),  Figure  3.8,  rad  vs.  time, 

which  has: 

(1)  Actuator  deflection  angle  0.05  rad. 
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(2)  Initially  and  within  the  0.28  sec  delay  time  a  small  transient  occurs 
with  an  insignificant  fluctuation. 

(3)  An  overshoot  1.3*  and  a  2.0  sec  settling  time. 

The  above  analysis  has  revealed  that  the  perfomance  of  the  system  is 
according  to  the  previous  mentioned  requirements. 

Also  should  be  mentioned  here  that  according  to  the  Eigenvalues  found 
executing  the  OPTSYS  program,  the  continuous  open  loop  system  Is  stable, 
since  the  s- plane  poles  are: 


s,  =  -159.724 

♦  j  18.9921 

s2  =  -159.724 

-  j  18.9921 

S3  =  -8.29048 

♦  j 8.05932 

S4  =  -8.29048 

-  j 8.05932 

S5  =  -3.75925 

♦  j  2.5 1463 

s§  =  -3.75925 

-  j  2.5 1463 

s7  =  -0.14393 

♦  JO 

S«  =  -0.00000425 

♦  jo 

S3  =  -0.04593 

♦  j 6.65959 

S|a=  -0.04593 

-  j 6.65959 

To  study  reconflgurable  flight  controls  the  missile  model  should 
be  capable  of  using  the  set  of  four  available  control  surfaces  to  produce  the 
required  forces  and  moments. 

To  achieve  this  each  control  surface  (  6tv  S2f  63,  S4 )  Is  considered 
as  an  Independent  component. 


The  equations  are  then  coupled  so  that  a  component  of  control 
surfaces  used  as  a  single  or  pair,  will  produce  the  appropriate  moments.  For 
example  deflected  at  an  certain  angle,  while  the  other  three  remained  at 
a  zero  angle,  will  produce  primarily  a  rolling  moment,  and  deflected 
together  with  $2  no  roll  effect  is  present,  as  we  will  see  in  Chapter  v. 

Taking  under  consideration  the  convention  control  surfaces  theory 
of  Chapter  II.,  for  the  B.T.T.  missiles,  and  for  the  configuration  shown  in  the 

Figure  3.4,  one  can  easily  show,  for  the  continuous  pitch  autopilot  that  the 
input  command  in  the  actuotor  network  (6^)  can  be  split  up  Into  independent 

components,  using  the  following  equation  which  compute  the  deflection 
Inputs  to  Equation  (  IH.B.4-12  ).  (Also  refer  to  Figures  2.2  and  3.4  which 
show  the  control  surface  positions  and  corresponding  positive  deflections). 

♦6pc  =  (-6,  +«4 )  f  2  (III.B.5- 1 ) 

Replacing  the  Input  vector  In  Equations  (III.B.4-1 1),  (III.B.4-12) 
and  (III.B.4-1 3)  with  the  above  Equation  and  applying  the  superposition 
principle,  gives  the  following  modified  airframe  equations: 

1 

6p  =  - 1 80.46p  ♦  168.4  —  (-  5,  ♦  82  -  63  ♦  64)  (III.B.5-2) 

2 

1  1 

—  (*  ♦  ^2  "  *  £4)  =  ~0.143  —  (-6^  ♦  ~  63  ♦  64)  ~  2.2308X2 

2  2 

♦  2.2308X4  -  340.58x1 0-*  X?  ♦  460.839x10-3^, 

-  460.839x1 0-3Xt  ♦  340.58x10-3X5  -  20.1845p 

(III.B.5-3) 


rfco  =  -  (X7  -  887. 1 24)6p  ♦  886.422 — (-6,  *  S2  ~  63  ♦  54  )  (III.B.5-4) 

2 


In  Figure  3.6  the  modification  of  the  commanded  tail  incidence  is 


presented. 
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Figure  3.6  (a)  Original  System;  (b)  Modified  System. 


Rearranging  the  above  equations  and  using  the  transfer  functions 
of  the  system  for  the  uncoupled  pitch  channel.  Equations  (III.B.4-4  to 
III.B.4-13),  the  following  equations  are  derived  and  can  be  modeled  in  a 
thirteenth  order  system  for  the  same  channel: 


X,  =  -150X,  ♦  150nz 

X2  =  X3  -  1.3531nZc+  1.353 IX, 

X3  =  -5X3  ♦  6.572nZc  -  6.572X, 


(III.B.5-5) 


(III.B.5-6) 


(III.B.5-7) 


*4  =  X5  *  59.26376p  (III.B.5-8) 

X5  =  - 1 49. 1 86x  1 0‘3X5  -  44.3557X4  ♦  1 .78875p  (III.B.5-9) 

Xs  =  X?-701.92x10-36p  (III.8.5- 1 0) 

X?  =  - 1 49. 1 86X 1 0-3X?  -  44.3557X*-  1 043.366?  (III.B.5- 1 1 ) 

5?  =  -188.46?  -  94.26,  ♦  94.252  -  94.25,  ♦  94.264  (III.B.5- 12) 


l,  =  -0.1436,  ♦  1.1 154X2  -  1.1 154X4  ♦  0.17029X3 

-0.23042^  ♦  0.23042X,  -  0.17029X5  ♦  10.0296p  (III.B.5- 13) 

&2  =  “0. 1 4362  -  1.1 154X2  ♦  1.1 154X4  “  0.17029X3 

♦  0.23042nzc  -  0.23042X2  ♦  0.17029XS  -  10.0926p  (III.B.5- 14) 

£3  =  -0.1436,  ♦  1.1 154X2  -  1.1 154X4  ♦  0.17029X, 

-  0.23042n2c  ♦  0.23042X,  -  0.17029X5  ♦  10.0296p  (III.B.5- 15) 

=  -0.1 4364  -  1.1 154X2  ♦  1.1 154X4  *  0.17029X3 

♦  0.2304211a.  -  0.23042X,  ♦  0. 1 7029X5  -  1 0.0296p  (III.B.5- 1 6) 

1^  =  -X7  ♦  887.1246p  ♦  443.21 16,  -443.21 162  ♦  443.21 16, 

-443.21164  (III.B.5- 17) 

We  thus  have  a  multiple-input-multiple-output  (MIMO)  system 
with  the  following  order 

•  X  :  1 3x  1 ,  the  time  derivative  of  X  matrix. 

•  F  :  13x13,  the  open  loop  dynamics  matrix  or  plant. 

•  X  :  13x1,  the  variables  asslgnements  matrix. 

•  G  :  13x5,  the  control  distribution  matrix. 

•  U  :  5x1,  the  control  Input  matrix. 

The  above  plant  control  and  output  matrices  are  presented  at  the 
Appendix  A  of  this  thesis. 
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Applying  the  same  step  Input  function  as  In  section  III.B.4,  which 
represents  the  a1  gee  command*,  and  zero  degrees  deflection  command  for 
the  four  control  surfaces  (1.e.,  61c,  S2c,  63c  64c  =  0)  the  time  response  plots 

of  the  following  state-varlahles  of  Interest  are  obtained  and  presented  in 
the  Figures  3.10  through  3.13. 

The  performance  of  the  system  (  q,  6p )  as  shown  In  Figures 
3.10  to  3.12  is  the  same  as  that  shown  in  Figures  3.6  through  3.8. 

In  Figure  3.13  Is  shown  the  time  response  of  the  control  surfaces 
deflection  angle,  for  the  Input  of  an  acceleration  equal  to  1  gee  command. 

Deflections  (  6U  S2,  h*  $4  °r  D's  ),  Figure  3.13,  rad  vs.  time, 
which  has: 

(1)  The  S,  and  63  have  the  same  sign  (negative),  as  it  was  expected 
according  to  the  sign  convention  of  Figure  2.2. 

(2)  The  S2  and  S*  have  the  same  response  as  above  but  with  a  positive 
sign. 

(3)  Insignificant  overshoot. 

(4)  A  2  seconds  settling  time. 

In  the  same  Figure  3.13  the  deflection  of  the  control  surfaces  Is 
Illustrated.  The  presented  deflection  Is  the  required  position  of  each 
surface  in  order  to  achive  the  commanded  one  gee  acceleration  (downwards). 

We  procceed  now  to  the  analysis  of  the  impaired  control  surfaces 
for  the  modified  system. 
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Figure  3.9  Pitch  Tall  Incidence  vs.  Time;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  o,=0;  Step  Input  =  1  g's. 
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Figure  3.10  Pitch  Normal  Acceleration  vs.  Time;  Uncoupled  Pitch  Channel; 

Unimpaired  Open  Loop  System;  a*=0;  Step  Input:  n^sl  g's  and 
*10  *50  *4C  =  0  degrees. 
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Figure  3.1 1  Pitch  Angular  Rate  vs.TIme;  Uncoupled  Pitch  Channel; 

Unimpaired  Open  Loop  System;  a,=0;  Step  Input:  n^sl  g's  and 

Sic.  $2*  $3c  $4c  =  0  degrees. 
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Figure  3.12  Pitch  Toll  Incidence  vs.  Time;  Uncoupled  Pitch  Channel, 

Unimpaired  Open  Loop  System;  %=();  Step  Input:  n,c=  1  g's  and 

5ic»  ^2c»  ^3c»  64c  =  0  degrees. 


C.  CONTROL  EFFECTIVENESS  DUE  TO  THE  CONTROL  SURFACES 

IMPAIRMENT  AND  COMPARISON 

The  time  response  plots  of  the  system  are  obtained  In  this  section 
for  two  types  of  damages.  The  first  type  Is  the  locked  or  hard  over 
surface  at  a  small  angle  and  the  second  type  Includes  the  Inoperatlvlty  of 
one,  two  or  three  actuators.  The  second  type  of  damage  could  also 
represent  the  destruction  of  the  control  surface(s),  i.e.v  shot  away. 

Note  that  the  first  phase  of  the  design  approach  for  the  CBTT 
autopilot  was  to  design  each  channel  independently  with  all  coupling 
between  channels  removed.  Sufficient  high  frequency  attennuatlon  was 
added  for  actuator  and  missile  design  so  that  the  resulting  missile  body 
angular  rates  and  control  surfaces  motion  would  represent  a  practical 
missile  design.  A  relationship  was  established  among  the  relative  speeds 
of  response  of  the  uncoupled  channels  In  order  to  meet  CBTT  requirements. 
The  cross-coupling  due  to  control  surfaces  Impairment  between  the 
longitudinal  and  lateral  motion  of  the  missile  Is  Ignored  when  each  control 
channel  Is  treated  Independently.  Cross-coupling  of  the  controls  is  thus 
considered  in  Chapter  V. 


One  control  surface  (6t)»  was  selected  to  be  deflected  or  hard 
over,  l.e.,  It  remained  at  a  fixed  angle  while  the  other  three  control 
surfaces  have  a  variable  deflection  angle. 

The  control  Input  was  a  step  Input  of  one  gee  In  the  normal 
acceleration  as  in  the  original  and  unimpaired  models.  For  the  damaged 
control  surface  (6t)  a  hard  over  deflection  of  0.2  rad  was  considered.  The 


time  response  plots  of  the  state-variables  of  Interest  were  obtained  and 
analyzed  with  the  performance  criteria  as  defined  in  section  (III.B.4): 

a.  Pitch  normal  acceleration  (nz  or  Nz),  Figure  3.14,  g's  vs. 
time,  which  has: 

(1)  Initially  and  with  a  fast  response,  the  missile  body  tends  to  go 
upwards  reaching  a  0.6  gee  acceleration.  This  Is  caused  by  the  fact 
that  at  the  beginning  of  the  command  the  damaged  surface  was 
deflected  at  a  +0.2  rad  angle  of  deflection  while  the  other  three 
were  at  their  steady  position.  Then  the  missile  goes  downwards 
until  the  -1.25  gee  acceleration  is  achieved.  This  Is  due  to  the 
deflection  of  the  unimpaired  surfaces  which  are  commaned  In  such  a 
way  as  to  counteract  the  damaged  (5t)  surface.  From  this  point  (0.2 
sec)  the  acceleration  rises  constantly  In  order  to  achieve  the  one 
gee  command. 

(2)  A  rise  time  0.75  seconds  and  a  negligible  overshoot. 

(3)  A  delay  time  about  0.5  sec 

(4)  Steady-state  error  0.005%  which  is  Insignificant. 

b.  The  pitch  angular  rate  (q  or  Q),  Figure  3.15,  rad/sec  vs.  time, 
which  has: 

(1)  Initially  and  due  to  the  fact  that  Is  locked  at  a  positive  angle  of 
0.2  rad,  to  achieve  the  acceleration  response,  the  required  body 
angular  rate  the  first  tenth  of  the  second  increases  very  fast  to 
0.40  rad/sec.  Then  It  begins  to  decrease  crossing  the  zero  rate  at 
0.2  sec,  corresponding  to  the  max  deceleration  (-1.25  gee's).  Then 
it  changes  direction  for  the  next  0.2  sec  until  the  value  of  -0.30 
rad/sec  Is  achieved.  From  this  point  the  rate  Is  Increased  again 
until  the  constant  value  of  -0.01  rad/sec  Is  achieved  at  the 
settling  time  1.4  sec. 

(2)  The  initial  increase  of  the  rate  is  due  to  the  produced  pitch  moment 
caused  by  the  deflection  of  the  damaged  control  surface  while  the 
other  three  are  at  their  steady  position. 


c.  Pitch  tall  incidence  (6p  or  DP),  Figure  3.16,  rad  vs.  time 
which  has: 

(1)  To  achieve  the  acceleration  response  and  due  to  the  fact  that  is 
locked  at  an  angle  of  0.2  rad,  initally  more  control  surface 
deflection  is  required  (  -0.13  rad).  At  0.2  sec,  when  the 
acceleration  begins  to  increase  and  the  angular  rate  is  zero,  the  tail 
incidence  has  the  positive  angle  of  0.03  rad  and  a  fluctuation  begins 
before  being  stabilized  at  a  constant  angle  of  0.06  rad. 

(2)  Settling  time  equal  to  1.8  sec. 

d.  Deflections  (  S„  S2,  S3,  S4  or  D's  ),  Figure  3.17,  rad  vs.  time. 

This  Figure  illustrates  the  time  response  of  the  deflection 

of  the  four  control  surfaces, when  one  of  them  (St)  is  locked  at  an  angle  of 
0.2  rad,  in  order  to  achieve  the  acceleration  response.  Initially  while  S2, 
S3  and  S4  control  surfaces  have  zero  angle  of  deflection,  S1  is  already  at 
the  0.2  angle.  Then  S2  and  S4  are  deflected  at  a  positive  angle  0.17  rad  and 
S3  at  a  negative  0.17  rad,  according  to  the  sign  convention  shown  in  Figure 
2.2. 

A  comparison  between  Figure  3.17  and  3.13  of  the  unimpaired 
system,  shows  that  to  achieve  the  acceleration  response  the  impaired 
system  requires  more  control  surface  deflection. 

Finally  to  cover  completely  this  class  of  impairment,  the  case  of 
a  negative  angle  (-0.2  rad)  at  the  St  surface  is  examined.  The  time 
response  of  the  same  state  variables  are  obtained  and  analyzed. 

a.  Pitch  normal  acceleration  (nz  or  Nz),  Figure  3.18,  g's  vs. 
time,  which  has: 

(1)  Initially  and  with  a  fast  response,  the  missile  body  tends  to  go 
downwards  reaching  a  0.6  gee  acceleration.  This  is  caused  by  the 


fact  that  at  the  beginning  of  the  command  the  damaged  surface  was 
deflected  at  a  -0.2  rad  angle  or  deflection  while  the  other  three 
were  at  their  steady  position.  Then  the  missile  goes  upwards  until 
the  4-1.75  gee  acceleration  Is  achieved.  The  produced  overshoot  is 
due  to  the  deflection  of  the  unimpaired  surfaces  which  are 
commanded  in  such  a  way  as  to  counteract  the  damaged  (6t)  surface. 
From  this  point  (0.3  sec)  the  acceleration  deacreases  constantly  in 
order  to  achieve  the  one  gee  command. 

(2)  A  small  0.18  second  rise  time  (reason  for  the  big  overshoot). 

(3)  A  delay  time  about  0.1  sec 

(4)  Steady-state  error  0.005%. 

b.  The  pitch  angular  rate  (q  or  Q),  Figure  3.19v  rad/sec  vs.  time, 
which  has: 

(1)  Initially  and  due  to  the  fact  that  5,  is  locked  at  a  negative  angle  of 
-0.2  rad  the  required  body  angular  rate  the  first  tenth  of  the  second 
decreases  very  fast  to  0.68  rad/sec.  Then  It  begins  to  increase 
crossing  the  zero  rate  at  0.3  sec,  corresponding  to  the  max 
deceleration  (1.75  gee's),  and  ft  changes  direction  for  the  next  0.5 
sec  until  the  value  of  0.13  rad/sec  is  achieved.  From  this  point  the 
rate  Is  decreased  again  until  the  constant  value  of  -0.01  rad/sec  is 
achieved  at  1.8  sec  settling  time  . 

(2)  The  big  overshoot  of  the  rate  Is  due  to  the  produced  pitch  moment 
caused  by  the  deflection  of  the  damaged  control  surface  while  the 
other  three  are  at  their  steady  position. 

c.  Pitch  tall  Incidence  (6p  or  DP),  Figure  3.20,  rad  vs.  time 
which  has: 

(1)  To  achieve  the  acceleration  response  and  due  to  the  fact  that  St  is 
locked  at  a  negative  angle  of  -0.2  rad,  Inltally  more  control  surface 
deflection  Is  required  (about  0.16  rad).  At  0.3  sec,  when  the 
acceleration  begins  to  decrease  and  the  angular  rate  Is  zero,  the  tall 
Incidence  has  the  positive  angle  of  0.01  rad.  A  fluctuation  begins 
before  the  deflection  being  stabilized  at  a  constant  angle  of  0.06 
rad. 
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(2)  Settling  time  equal  to  1.8  sec. 

d.  Deflections  (  8U  S2,  S3,  84  or  D's  ),  Figure  3.21,  rad  vs.  time. 

This  Figure  Illustrates  the  time  response  of  the  deflection 
of  the  four  control  surfaces,  when  one  of  them  (6))  Is  locked  at  a  negative 
angle  of  -0.2  rad.  In  order  to  achieve  the  acceleration  response.  Initially 
and  when  the  S2,  8?  and  64  control  surfaces  have  zero  angle  of  deflection 
St  Is  already  at  the  -0.2  angle.  Then  S2  and  84  are  deflected  at  a  negative 
-0.09  rad  angle  and  83  at  a  positive  0.09  rad. 

A  comparison  between  Figure  3.21  and  3.17  indicates  that  when 
the  surface  S,  is  hard  over  at  a  negative  angle: 


(1)  The  undamaged  control  surfaces  In  order  to  counteract  the  failure 
and  achieve  the  commanded  acceleration,  are  deflected  at  a  smaller 
angle. 


(2)  The  undamaged  control  surfaces  are  commanded  to  turn  at  angles 
with  opposite  sign  comparing  with  those  in  Figure  3.17  (S, 
positively  hard  over). 


The  5t  control  surface  was  selected  to  simulate  the  inoperative 
actuator.  The  corresponding  state  variable  S,  was  therefore  set  equal  to 
zero.  The  step  input  of  one  gee  remained  the  same.  The  Initial  condition  on 
the  control  surface  was  set  at  zero.  This  type  of  damage  could  also 
represent  the  distinction  of  the  control  surface. 

For  the  above  case  the  time  response  plots  of  the  state  variables 
of  interest  were  obtained  and  analyzed  according  to  the  defined 
performanse  criteria  In  Sectlor  (III.B.4).  Later,  In  Section  5,  a  comparison 
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Figure  3.15  Pitch  Angular  Rate  vs.Time;  Hard  Over  6t  Control  Surface; 
at=0;  Step  Input:  n^sl  gee,  81c=0.2  rad,  62»  ho  $4c=0- 
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Figure  3.18  Pitch  Normal  Acceleration  vs.  Time;  Hard  Over  Sy  Control 
Surface;  a,sO;  Step  Input:  n^l  gee,  S1c=-0.2  rad,  $2, 
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Figure  3.19  Pitch  Angular  Rate  vs.TIme;  Herd  Over  8,  Control  Surface; 

a,=0;  Step  Input:  n^l  gee,  Sur-0-2  rad,  82c,  5JC,  8^=0. 
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between  the  different  types  of  damages  and  the  unimpaired  system  was 
performed. 

.  a.  Pitch  normal  acceleration  (rig  or  Nz),  Figure  3.22f  g's  vs. 
time,  which  has: 

(1)  Time  constant  or  rise  time  0.41  seconds. 

(2)  A  0.005%  steady  state  error. 

(3)  Minor  overshoot  and  delay  time  0.35  sec  and 

(4)  A  1.6  sec  settling  time  at  which  the  acceleration  response  is 
achieved. 

b.  Pitch  angular  rate  (q  or  Q),  Figure  3.23,  rad/sec  vs.  time, 
which  has: 

(DA  maximum  body  angular  rate  of  0.1?  rad/sec  Is  achieved  at  0.3 
second. 

(2)  Minor  overshoot  and  a  settling  time  equal  to  1.6  after  which  the 
angular  rate  remains  constant  at  -0.01  rad/sec.  This  Is  due  to  the 
fact  that  the  transfer  function  q/8p  (Equation  III.B.3-1)  Is  a  type 
zero  system  and  has  a  finite  error  to  a  step  Input. 

c.  Pitch  tall  Incidence  (6p  or  DP),  Figure  3.24,  rad  vs.  time, 
which  has: 

(1) A  0.054  rad  command  from  the  actuator  networK. 

(2)  A  0.4  sec  delay  time  and  negligible  overshoot. 

d.  Deflections  (St,  S2,  S3,  S4  or  D's),  Figure  3.25,  rad  vs.  time, 
which  has: 

(DA  positive  deflection  angle  of  0.035  rad  for  the  S2  and  $4  actuators 
while  the  S3  Is  deflected  at  a  negative  and  equal  angle 
(-0.035  rad),  and  8)  as  damaged  remains  constantly  at  zero. 
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(2)  A  1.6  sec  settling  time  for  each  of  the  unimpaired  surfaces. 


(3)  A  rise  time  equal  to  0.7  sec  and  minor  overshoot  correspondingly 
for  the  undmaged  surfaces. 

(4)  The  deflection  angle  of  the  unimpaired  control  surfaces  has  a 
pitching  moment  effect  on  the  system. 

A  comparison  between  the  deflection  angle  In  Figures  3.13,  3.17, 
3.21,  and  3.25  Indicates  that  for  the  same  step  Input: 

(1)  When  one  control  surface  Is  Inoperative  or  shot  away  the  system 
commands  the  undamaged  surfaces  at  larger  angles  than  those 
obtolned  with  the  unimpaired  system. 

(2)  When  one  control  surface  Is  Inoperative  or  shot  away  the  system 
commands  the  undamaged  surfaces  at  smaller  angles  (In  comparison 
with  the  hard  over  type  of  damage). 

3.  Damage  Two  Control  Surfaces 

In  this  case  two  control  surfaces,  on  the  same  side  of  the 
missile  were  considered  damaged  (St  and  $2).  The  two  state  variables  S, 
and  S2  were  therefore  set  equal  to  zero.  The  step  input  of  one  gee 
remained  the  same.  The  initial  condition  on  the  control  surfaces  was  set 
at  zero.  Then  the  time  response  plots  of  the  state  variables  were  obtained 
and  analyzed. 

a.  Pitch  normal  acceleration,  Figure  3.26,  which  has  a  0.6  rise 
time,  0.38  sec  delay  time  and  0.005X  steady  state  error.  Note  the  faulty 
behavior  at  the  start  of  the  response  (0  to  0.15  sec).  Because  the  system 
is  a  nonminimum  phase  system  (i.e.,  a  zero  lies  in  the  right-half  s  plane) 
the  transient  response  starts  out  In  the  opposite  direction  to  the  input 
i.e.,  there  is  an  inherent  delay  in  the  transient  response  of  the  system. 


b.  Pitch  angular  rate,  Figure  3.27,  which  has  a  maximum  rate 
of  -0.16  rad/sec  at  0.3  sec.  Then,  increases  until  the  constant  negative 
rate  (-0.01  rad/sec)  Is  achieved  at  1.3  seconds. 

c.  Pitch  tall  Incidence,  Figure  3.26,  which  has  a  0.8  sec  rise 
time  and  0.3  sec  delay  time.  The  command  from  the  actuator  network  is 
0.054  rad  which  is  achieved  at  1.6  sec  settling  time. 

d.  Deflections  (6t>  S2,  6 3,  S4  or  D's),  Figure  3.29,  In  which 
the  two  damaged  control  surfaces  (St  and  S2)  remain  always  at  zero 
(damaged)  and  the  two  others  opposite  deflected  at  the  angle  of  0.053  rad 
angle  after  a  lapse  time  of  1.6  sec.  Note  that  for  this  system  with  two 
Impaired  control  surfaces  more  angle  of  deflection  Is  required  (In 
comparison  with  the  system  having  one  actuator  shot  away). 

4.  Damage  Three  Control  Surfaces 

In  the  last  case  three  of  the  control  surfaces  Su  S2  and  S3  were 
damaged  or  inoperative. 

With  a  step  Input  of  one  gee's  for  the  commanded  normal 
acceleration  as  In  the  unimpaired  and  other  type  damaged  models,  the 
plots  of  the  state  variables  of  Interest  were  obtained  and  analyzed. 

Figures  3.30  through  3.33  illustrate  the  response  of  the  output 
variables.  Figure  3.30  shows  that  the  commanded  one  gee  acceleration  can 
still  be  achieved  In  the  absense  of  three  control  surfaces  (61t  s2,  S3). 
However,  some  oscilatlon  Is  evident  Into  the  response.  The  loss  of  the 
surfaces  has  produced  a  more  lightly  damped  system  as  evidenced  in 
Figures  3.30  through  3.33.  The  oscillatory  and  convergent  behavior  of  the 
pitch  rate  and  the  undamaged  surface  (84),  shown  In  Figures  3.31  and  3.33, 


Illustrate  how  the  surface  attempt  to  'fill  in'  for  the  absence  of  the 
Inoperative  control  surfaces.  Specifically  the  response  of  the  pitch  normal 
acceleration  (Figure  3.30)  has: 


(1)  A  1.01  sec  time  constant. 

(2)  A  0.38  sec  delay  time  and 

(3)  A  3.0  sec  settling  time  at  which  the  steady-state  value  of  the 
commanded  acceleration  is  achieved. 

A  comparison  with  the  response  of  the  corresponding  deflection 
of  the  unimpaired  system  (Figure  3.13)  and  the  other  inoperative  cases. 
Figures  3.21,  3.25  and  3.29,  Indicates  that  as  the  number  of  the 
Inoperative  control  surfaces  Increases,  the  deflection  of  the  undamaged 
actuatorfs)  also  increases  and  more  time  Is  required  to  achieve  the 
steady-state  value. 


In  this  section  a  performance  comparison  of  each  type  of 
Impairment  with  the  undamaged  system  is  performed  and  an  analysis  of 
their  corresponding  time  response  plots  is  presented. 

In  order  to  use  the  N.P.S.  EA5YPL0T  program  the  state  variables 
had  to  be  input  as  data.  This  was  accomplished  by  writing  a  computer 
program  ( Appendix  0). 

With  the  above  data  Figure  3.34  is  obtained.  This  corresponds  to 
pitch  normal  acceleration  (n,  or  NZ)  for  the  unimpaired,  hard  over  Slf  and 

inoperative  or  shoot  away  one  or  more  control  surfaces,  configurations. 

The  hard  over  case  shows  greater  excursion  In  the  normal 
acceleration  with  a  peak  value  of  -1.25  gee's  at  0.2  sec  being  achieved. 
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□  NZ, PITCH  NORMAL  ACCELERATION-G'S 
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Figure  3.22  Pitch  Normal  Acceleration  vs.  Time; 

Inoperative  Sy  Control  Surface;  a^O; 

Step  Input:  n2c  =1  gee,  S1c  =  0.2,  82c,  S3c,  s«c  =  0  rad. 
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Figure  3.24  Pitch  Tall  Incidence  vs  Time; 

Inoperative  Control  Surface;  a^O; 

Step  Input:  n2c  =1  gee,  S1c  =  0.2,  S2»  ho  s4c  =  0 
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Figure  3.27  Pitch  Angular  Rate  vs  Time; 

Inoperative  Sj  and  S2  Control  Surfaces;  a^O; 

Step  Input:  n2c  =1  gee,  $1c  =  0.2,  $2c,  5^  =  0  rad. 
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Figure  3.28  Pitch  Tell  Incidence  vs  Time; 

Inoperative  and  S2  Control  Surfaces;  a^O; 

Step  Input:  nZc  =1  gee,  5,c  =  0.2,  S2c>  53o  S4c  =  0  rad. 
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Figure  3.30  Pitch  Normal  Acceleration  vs.  Time; 

Inoperative  6U  S2  and  63  Control  Surfaces;  a^O; 

Step  input:  n2c  =1  gee,  51c  =  0.2,  62c>  he*  ^4c  =  0  rad. 


Q, pitch  angular  rate-rad/sec 
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Figure  3.31  Pitch  Angular  Rate  vs  Time; 

Inoperative  62  and  63  Control  Surfaces;  8e=0; 

Step  input:  nZc  =1  gee,  S1c  =  0.2,  S2c,  63o  5^  =  0  rad. 
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Figure  3.33  Control  Surfaces  Deflections  vs.  Time; 

Inoperative  61(  $2  and  53  Control  Surfaces;  a8=0; 

Step  input:  n2c  =1  gee,  S1c  =  >0.2,  S:r>  $3c,  s*  =  0  rad. 


The  Inoperative  cases  are  all  quite  simitar  except  for  the  last 
case  where  only  one  control  surface  Is  available.  In  this  latter  case  an 
oscillatory  behavior  Is  evident.  The  performance  criteria  to  the  one  gee 
command  are  summarized  In  Table  3.2. 

TABLE  3.2 


Sustem 

Unimpaired-  Hard  Over-  Inoperative  Control  Surfaces 

One  Two  Three 

Time  Constant(sec):  0.50 _ 0.4 1 . £>.60 _ .0.60... _ 1 .0 1 

Delay  Tlme(sec):  0.28 _ ..0.5 _  .J0.35......J0.38.........  0.40 

Overshoot(X):  2.2X _ jO.O  1 X _ 0.0 1  X.....0.0 1 X _ 0.0 1 X 

Figure  3.35  compares  the  unimpaired  system  with  the  three 
cases  of  control  surfaces  Impairments.  The  curves  include  the  following 
four  cases: 

(1)  Original  system 

(2)  Inoperative  one  control  surface 

(3)  Inoperative  two  control  surfaces 

(4)  Inoperative  three  control  surfaces 

This  Figure  shows  that  as  the  number  of  the  damaged  control 
surfaces  Increases  the  delay  time  Increases  also  from  0.28  sec  for  the 
unimpaired  system  to  0.41  sec  for  the  system  with  three  Inoperative 
control  surfaces.  The  Increase  also  of  the  time  constant  from  0.5  sec 


(unimpaired  system)  to  about  1  sec  (three  Inoperative  actuators), 
Indicates  that  as  the  number  of  the  damaged  control  surfaces  increases 
more  time  Is  required  to  achieve  the  steady-state  value  of  the  normal 
acceleration.  Also  observe  the  initially  faulty  response  of  the  nonminimum 
phase  system  to  the  step  Input  which  decreases  as  the  number  of  the 
inoperative  control  surfaces  Increases. 

Figure  3.36,  shows  that  to  achieve  the  acceleration  response  the 
system  with  the  hard  over  5j  type  of  damage,  requires  more  angular  rate 
(q  or  Q),  compared  with  the  unimpaired  and  the  other  three  damaged 
system  cases.  The  produced  pitching  moment  Is  evident.  Initially  and 
because  the  51  control  surface  Is  deflected  at  the  0.2  rad  angle  while  the 
other  three  are  at  the  zero  position,  the  first  tenth  of  the  second  the  rate 
Increases  from  zero  to  0.45  rad/sec.  The  rate  then  begins  to  decrease 
very  rapidity  to  a  rate  of  -0.38  rad/sec  at  0.4  sec. 

Figure  3.37  compares  the  unimpaired  system  with  three  cases  of 
control  surfaces  Impairments.  The  curves  Include  the  following  four  cases: 

(1)  Original  system 

(2)  Inoperative  one  control  surface 

(3)  Inoperative  two  control  surfaces 

(4)  Inoperative  three  control  surfaces 

This  Figure  shows  that  as  the  number  of  the  damaged  control 
surfaces  Increases,  the  angular  rate  of  change  Is  reduced,  and  more 
pronounced  fluctuations  occur,  1.e.,  the  system  becomes  more  oscillatory. 
These  fluctuations  are  connected  with  the  pitching  moment  induced  in  the 
missile  by  the  failed  actuators. 


Finally  Figures  3.38  and  3.39  compare  the  unimpaired  system 
with  four  cases  of  control  surfaces  impairments.  The  curves  include  the 
following  five  cases: 

(1)  Original  system 

(2)  Hard  over  one  control  surface 

(3)  Inoperative  one  control  surface 

(4)  Inoperative  two  control  surfaces 

(5)  Inoperative  three  control  surfaces 

Figure  3.38  and  3.39  show  that  the  required  (steady-state)  pitch 
tall  Incidence  for  the  uncoupled  longitudinal  autopilot  is  0.054  rad 
(unimpaired  or  damaged). 
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Figure  3.34  Comparison  of  the  Pitch  Normal  Acceleration  between  the 
Unimpaired  System  and  the  Four  Cases  of  Impairment. 
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Figure  3.36  Comparison  of  the  Pitch  Angular  Rate  between  the 

Unimpaired  System  and  the  Four  Cases  of  Impairment. 
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Figure  3.37  Pitch  Angular  Rate;  Uncoupled  open  loop  Channel: 
Unimpaired  and  Damaged  Control  Surfaces. 
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Figure  3.38  Comparison  of  the  Pitch  Tail  Incidence  between  the 
Unimpaired  System  and  the  Four  Cases  of  Impairment. 
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Figure  3.39  Pitch  Tall  Incidence;  Uncoupled  open  loop  Channel: 
Unimpaired  and  Damaged  Control  Surfaces. 
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In  this  chapter  the  design  and  analysis  of  the  continuous  uncoupled 
lateral  (1.e.,  yaw  and  roll)  channels  of  a  B.T.T.  autopilot  derived  In  [Ref.  1 
and  4]  is  reviewed.  A  modified  lateral  autopilot  Is  then  designed  with  the 
four  control  surfaces  split  Into  independent  components.  The  model  of 
elliptical  (yaw  channel)  and  circular  (roll  channels)  configurations  are  both 
shown  in  Figure  4.1.  Both  airframes  have  the  same  cross  sectional  area 
distribution  and  specifically,  the  circular  cross  sectional  body  has  a  closure 
ratio  Abase/ Amax  of  0.69  with  Amax  occurlng  at  68X  missile  body,  whereas 
the  elliptical  airframe  has  a  3:1  cross  section. 

A  comparison  between  the  two  systems  (original  and  modified)  shows 
that  they  are  Identical.  The  transient  response  plots  and  their  analysis 
demonstrate  the  performance  of  the  original  and  modified  yaw  and  roll 
channels. 

Having  established  the  validity  of  the  modified  system  one  then  uses 
the  multlple-input/multlple-output  control  theory  and  state-space 
representation  [Ref.3]  to  Investigate  the  performance  of  the  modified 
lateral  autopilot  for  various  classes  of  damages. 

For  the  Investigation  of  the  lateral  channel,  two  classes  of 
impairments  at  the  control  surfaces  are  considered: 

(1)  Locked  or  hard  over  and 

(2)  Inoperative  with  the  control  surfaces  centered  or  shot  away. 
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The  plots  of  the  unimpaired  modified  yaw  and  roll  autopilots  are  used 
as  a  basis  for  comparison.  The  impairment  induced  effects  on  roll  angle,  X 

rate  and  tall  incidence  are  assessed.  6 


B.  UNCOUPLED  LINEAR  YAW  CHANNEL  (ORIGINAL) 

1.  Aerodynamic  Model.  Airframe  Conflouratlon  and  Control  Law 

A  general  block  diagram  of  a  B.T.T.  autopilot  with  all  Its  channels 
is  shown  in  Figure  3.2.  Inertial  acceleration  commands  are  applied  In  polar 


coordinates,  i.e.,  magnitude  of  the  command  (rt)  applied  to  pitch  autopilot  ■ 

and  the  direction  (fc)  Is  applied  to  the  roll  autopilot.  The  yaw  autopilot  is 

slaved  to  the  roll  autopilot  In  order  to  minimize  sideslip  angle  by  y 

ii 

coordinating  the  missile  yaw  and  roll  motion.  Achieved  maneuver  plane  | 

accelerations  in  rectangular  coordinates,  (tfe)  and  (r^),  were  determined  by 

resolving  achieved  body-fixed  accelerations,  (r^)  and  (r^)  through  the  X 


missile  roll  angle  Op),  I.e.,  the  Euler  angles  6  and  ip  are  assumed  to  be  I 

sufficient  small,  so  that  the  sideslip  angle  (0)  is  minimized  and  the  motion  :: 

* 

between  the  yaw  and  roll  channels  is  coordinated. 

A  general  diagram  of  the  uncoupled  yaw  channel  where  the  I 

aerodynamic  and  yaw  control  law  are  Involved  Is  shown  In  Figure  4.2.  In 

Figure  4.3  [Ref.  1]  the  yaw  control  law  Is  presented.  The  yaw  normal  ^ 

acceleration  (i\yc)  Is  not  used  to  commond  the  coupled  B.T.T.  autopilot,  V 


although.  It  is  used  for  the  design  and  analysis  of  the  uncoupled  channel.  The 
command  used  by  the  coupled  system  is  shown  in  dashed  lines  and  is  a  yaw 
angular  rate  command  (rc).  The  yaw  control  tow  is  governed  by  missile  body 
yaw  ongulor  rote  (r)  and  yaw  normal  acceleration  ( ). 
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The  flight  conditions  are  the  same  as  was  mentioned  for  the  pitch 
channel,  (l.e.,  60  Kft  altitude  and  hack  number  3.95).  The  yaw  control  law 
determines  the  required  yaw  tall  Incidence  (6yc)  command  to  the  actuator 
network  which  Is  approximated  as  a  first-order  lag  at  30  Hz. 


2.  Requli 


lynuji 


Its  for  the  Classical  Design 


Requirements  for  the  classical  design  approach  are: 


a.  High  frequency  attenuation  in  actuator  command  branch  >  15  dB  at 
100  rad/sec  and  zert>  angles  of  attack  and  sideslip. 

b.  Relative  stability:  gain  margins  i  6  dB,  phase  margins  >30  deg  with 
a  goal  of  1 2  dB  and  50  degrees. 

c.  Acceleration  time  response: 

(1)  63%  time  constant  of  approximately  0.4  sec  for  a  step  command  of 
acceleration  (fy)  at  flight  conditions  of  interest  and  at  a,  and  0=0. 

(2)  Overshoot  i  10%  for  better  sideslip  control. 

(3)  Steady  state  error  not  necessarily  zero. 

3.  Transfer  Funciions  of  Aerodynamic  node)  and  Analyst? 

The  aerodynamic  transfer  functions  of  the  uncoupled  yaw  channel 
autopilot  are  according  to  Karadimas,  C.  thesis  [Ref.  4]: 
a.  Transfer  function  of  yaw  angular  rate: 

I 

0.15  { -  s  +  1} 

r  0.0748 

- = - (deg/sec/deg)  (IV.B.3- 1 ) 
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b.  Transfer  function  of  Yaw  Normal  Acceleration: 
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Rearranging  the  above  transfer  functions  In  terms  of  the  variable 
pairs  (r,£y)  and  (r^Sy)  respectively  and  applying  the  Inverse  Laplace 

transformation  a  set  of  linear  differential  equations  can  be  obtained. 
Then,  using  rules  of  state-space  representation,  the  following  equations 
are  obtained  [Ref.  4]: 


X,  =  X2  -  29.2321  Sy 

(IV.B.3-3) 

X2  =  -16.7224X,  -  0.0434X2  ♦  0.2687Sy 

(IV.B.3-4) 

2,  =  Z2  -  0.080  16y 

(IV.B3-5) 

t2  =  -16.7224!  -  0.0435Z2  -  150.77246Y 

(IV.B.3-6) 

Y  =  3.3574Z,  -4Y  ♦  6.434 1SY  -  3.3574%c 

$Yc=  6.06X,  ♦  0.606X2  ♦  2.04 13Z,  ♦  3.648Y  -13.8612 

(IV.B.3-7) 

-2.0413% 

(IV.B.3-8) 

$Y  =  188.4Syc  '  188.45y 

(IV.B.3-9) 

Using  state-space  representation,  the  Equations  (IV.B.3-3)  through 
(IV.B.3-9)  can  be  modeled  In  a  seventh-order  system  of  the  form: 
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with  the  continuous  plant  system  F  and  control  matrix  G  as  follows: 
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where  the  state  variables  are: 

-  X1(  X2  :  yaw  angular  rates 

-  Z, ,  Z2  :  yaw  normal  acceleration 

-  Y  :  output  of  acceleration  compensator 

-  $yc  :  Input  command  In  the  actuator 

-  Sy  :  yaw  tall  Incidence 


Using  a  step  Input  function  which  represents  the  “1  gee  commond' 
at  ae  =  0,  the  time  response  plots  (Figure  4.4),  the  Vaw  Normal 
Acceleration  (ity)  (Figure  4.5),  and  the  Vaw  Angular  Rate  (R)  and  the  Vaw 


Tall  Incidence  (6y)  (In  Figure  4.6)  were  obtained.  The  transient  response 

satisfies  the  performance  criteria  and  definitions  given  at  Section  (III.B.4), 
and  ore  Identical  with  the  results  presented  In  [Ref.  41. 

In  Figure  4.4  the  Yaw  Normal  Acceleration  Is  plotted  which  has: 

(1)  A  0.31  sec  time  constant. 

(2)  A  0.3  sec  delay  time  and  6.4%  overshoot. 

(3)  A  0.018  steady-state  error. 

(4)  Note  the  faulty  behavior  at  the  start  of  the  response  (0  to  0.2  sec). 
Because  the  system  Is  a  nominimum  phase  system  (l.e.,  a  zero  lies 
In  the  right-half  s  plane)  the  transient  response  starts  out  In  the 
oppslte  direction  to  the  Input  i.e.,  this  Is  a  Inherent  delay  In  the 
transient  response  of  the  system. 

In  Figure  4.5  the  Yaw  Angular  Rate  1e  presented  which  has: 

(1)  A  maximum  body  angular  rate  Of  0.6  rad/sec  achieved  at  0.3  second. 

(2)  A  minor  undershoot  and  2.2  sec  settling  time  after  which  the 
angular  rate  remains  constant  at  0.01  rad/sec.  This  Is  due  to  the 
fact  that  the  transfer  function  r/Sy  (Equation  IV.B.3-1)  is  a  type 
zero  system  and  has  a  finite  error  to  a  step  Input. 

Finally  Figure  4.6  shows  that  a  negative  angle  of  deflection  of 
-0.14  rad  is  required  In  order  to  achieve  the  commanded  step  input 
acceleration. 

The  above  analysis  has  revealed  that  the  performance  of  the  system 
confirms  to  the  previous  mentioned  requirements. 

From  the  pole-zero  plot  In  Figure  2.7  [Ref.4]  we  can  see  that  the 
continuous  open  loop  system  is  stable,  since  the  s-plane  poles  are: 

Sj  =  -175.010  -jO  (IV.B.3- 1 0) 
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0.3  1.0  1.3  2.0  2.3  3.0 
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Yaw  Tail  Incidence  vs.  Time;  Uncoupled  Yaw  Channel 
Autopilot;  Continuous  Open  Loop  System;  Elliptical  Airframe; 
Step  Input  HY  =  1  gee;  at  =  0. 


S2  =  -5.7665 

♦J  10.4804 

(IV.B.3- 1 1 ) 

Sj  =  -5.7665 

-J  10.4804 

(IV.B.3- 1 2) 

S4=  -0.0217473  +J4.08919 

(IV.B.3- 1 3) 

S5  =  -0.0217473 

-J4.08919 

(IV.B.3- 1 4) 

S6  =  -2.9484 

♦j 2.99489 

(IV.B.3-15) 

ST  =  -2.9484 

-j 2.99489 

(IV.B.3- 1 6) 

C.  UNCOUPLED  LINEAR  OPEN  LOOP  ROLL  CHANNEL  (ORIGINAL) 

1.  The  Aerodynamic  Model,  the  Airframe  Configuration  and  the 

Control  Law 

Generally,  the  purpose  of  a  roll  channel  autopilot  Is  to  command 
a  roll  Op)  signal.  In  order  to  roll  the  missile  through  an  angle  (9) 
measured  from  the  vertical  axis. 

While  the  pitch  channel  acceleration  Is  commanded  to  produce 
the  total  magnitude  of  the  guidance  acceleration  command,  the  roll 
channel  of  the  B.T.T.  autopilot  Is  commanded  to  roll  the  missile  so  as  to 
put  the  prefered  maneuver  direction  in  the  direction  of  the  guidance 
acceleration  command,  as  Is  shown  In  Figure  3.2. 

A  block  diagram  of  the  uncoupled  roll  channel  Is  shown  In  Figure 
4.7,  and  the  roll  control  law  In  Figure  4.8  (  both  Figures  taken  from 
[Ref.  H).  In  these  diagrams  the  same  aerodynamic  model  Is  used. 

As  is  shown  in  Figures  4.7  and  4.8  the  control  law  Is  commanded 
by  a  roll  angle  (fc)>  is  governed  by  roll  angular  rate  (P)  and  roil  angle  (<p) 
and  determines  the  required  command  (6*,)  to  an  actuator.  This  is 
approximated  as  a  first  order  lag  at  188.4  rad/sec. 


The  circular  airframe  configuration  Is  chosen  for  this  channel 
and  Is  shown  in  Figure  3.1.  The  flight  conditions  are  the  same  as  for  the 
Yaw  and  Pitch  channels,  1.e.  60000  ft  and  Mack  number  3.95. 


2.  Requirements  for  the  Classical  Design 

Requirements  for  the  classical  design  approach  are: 

a.  High  frequency  attenuation  in  Actuator  command  branch  i  15  dB  at 
100  rad/sec  and  zero  angle-of -attack.  This  requirement  limits  the 
speed  of  roll  angle  response. 

b.  Relative  stability:  gain  margins  i  30°  with  a  goal  of  12  dB  and  50°, 
i.e.  the  same  as  for  the  yaw  channel. 

c.  Acceleration  time  response : 

(1)  63*  time  constant  of  0.5  seconds  for  a  step  command  of  roll 
angle  at  the  flight  condition  of  interest  and  zero  angle  of  attack. 

(2)  Overshoot  i  10*.  This  requirement  is  conservative  and  may  be 
relaxed  depending  upon  guidance  results. 

(3)  Zero  steady  state  roll  angle  error.  If  not  zero,  it  will  affect 
guidance. 

3.  Transfer  Funclons  of  the  Aerodunamlc  Model  and  Analysis 
The  aerodynamic  transfer  function  of  the  uncoupled  roll  channel 

autopilot  are  [Ref.  4]. 

a.  Transfer  functions  of  the  Roll  Angular  rate: 
li  =  (q'Sd/Ixx)x  57.3  x  C|  jr  Sr  (rad/sec2) 

where: 

q*  =  1650  (lbs/ft2) 

3=3.  14159  ft2 


(IV.C.3- 1 ) 


Ixx  =40  slug-ft2 
C,  jD  =0.03 1 
£  =  8.0346  6r 

b.  Transfer  function  of  Roll  Angle: 

9  =  P  (rad/sec)  (IV.C.3-2) 

Rearranging  the  above  equations  along  with  those  given  in 
{Ref.  1]  for  the  actuator  and  using  rules  of  state-space  representation,  the 
following  equations  are  obtained  [Ref.  4]  and  verified: 


P  =  8.034628b  (IV.C.3-3) 

9  =P  (IV.C.3-4) 

X  = -17.69 -8X+1 7.69*  (IV.C.3-5) 

Y  =  -0.75325P  -0.88359  +0.35 165X  -5Y  -0.403345R 

♦0.883529*  (IV.C.3-6) 

X,  =  -6X,  ♦0.628229  18d  (IV.C.3-7) 

iRo=  -0.1 027 IP  -0.120489  +0.04795X  +14.3182Y 

- 1 4. 1 8 1 84X,  - 1  58b*  (IV.C.3-8) 

=  10795.328**  -188.48r  (IV.C.3-9) 


Utilizing  state-state  representation,  the  above  Equations  can  be 
modeled  as  a  seventh-order  system  in  the  format  X  =  FX  ♦GU  with  the 
continuous  plant  system  F  and  control  matrix  G  as  shown  in  Table  4. 1  and 
the  state  vector  with  the  corresponding  state  variables  as  follows. 
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P  Roll  angular  rate 

<P  Roll  angle 

X  Output  of  roll  angle  compensator 

Y  Output  of  roll  rate  compensator 

Xf  Ouput  of  pseudodlfferenstator 

fa  Input  command  tn  the  actuator 

Sr  Roll  tall  incidence 

Executing  the  Interactive  Implementation  of  the  optimal  systems 
control  design  program  (OPTSYSX)  using  an  Input  step  function  which 
represents  “1  radian  command’  at  zero  trim  angle  of  attack,  the  graphical 
time  response  plots  were  obtained  for  the  state  variables  of  Interest. 

Figure  4.9  shows  the  Roll  Angular  Rate  (P)  (rad/sec  vs.  time). 
Figure  4.10  the  Roll  Angle  (f)  (rad  vs.  time),  and  Figure  4.1 1  the  Roll  Tall 
Incidence  (Sr)  (deg.  vs.  time).  A  close  observation  In  Figure  4.10  Indicates 

that  for  the  commanded  Input  vector  the  time  response  plot  of  the  roll 
angle  has  a  0.55  sec  time  constant,  5.3X  overshoot  and  a  steady-state 
error  equal  to  zero.  These  results  are  In  accordance  with  the  requirements 
mentioned  In  the  previous  section. 

From  the  pole-zero  plot  of  Figure  2.23  (Ref.4)  we  can  mention 
here  that  the  continuous  open  loop  system  Is  stable,  since  the  s-plane 
poles  are: 


s,  =  -174.785  -jO 
s2  =  -9.2510  +j 28.4098 


(IV.C.3- 1 0) 
(IV.C.3- 1 1 ) 


8 


-9.25097  -j 28.4098 
-2.46608  ♦J2-71 152 
-2.46608  +J2.71 152 
-8.98209  *J0 
-5.19932  +J0 


(IV.C.3- 1 2) 
(IV.C.3- 1 3) 
(IV.C.3- 1 4) 
(IV.C.3-15) 
(IV.C.3-16) 


To  study  reconfigured  e  flight  controls  the  missile  model  should 
be  capable  of  using  the  set  of  four  available  control  surfaces  to  produce  the 
required  forces  and  moments. 

To  achieve  this  each  control  surface  (51t  S2.  h»  was  considered 
as  an  independent  component. 

The  equations  were  then  coupled  so  that  a  component  of  control 
surfaces  used  alone  or  in  conjunction  with  another  control  surface,  will 
produce  the  appropriate  moments. 

Talcing  under  consideration  the  conventional  control  surface  theory 
of  Chapter  II.  for  the  Bank-to-Turn  missiles,  and  for  the  circular 
configuration  shown  in  Figure  3.4,  one  can  easily  show,  for  the  continuous 
open  loop  roll  autopilot  that  the  input  command  in  the  actuator  network 
(5rc)  can  be  split  up  Into  Independent  components,  using  the  following 
equation  which  computes  the  deflection  Inputs  to  Equation  (IV.C.3-8).  Also 
refer  to  Figures  2.2  and  3.4  which  show  the  control  surface  position  and 
corresponding  positive  deflections. 


♦6rc  =  "  $2  “  $3  *$4)  f  4 


(IV.C.4- 1 ) 
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Replacing  the  Input  vector  in  Equation  (IV.C.3-8)  and  (IV.C.3-9) 
with  the  above  Equation  and  applying  the  superposition  principle,  gives  the 
following  modified  airframe  equations: 

=  -188.46g  ♦  10795.32  (6t  -  62  -  63  ♦  64)  /  4  (IV.C.4-2) 

(4,  -$2  -4,  /  4  =  -0.1 027 IP  -0.12048*  +0.04795X  +14.318Y 

-1418 IBXi  -15/4  -  0.140676r 

-  0.12048*c  (IV.C.4-3) 

In  Figure  4.6  the  modification  of  the  commanded  tail  incidence  is 
presented. 


Figure  4.6  (a)  Original  System;  (b)  Modified  System. 

Rearranging  the  above  equations  and  using  the  transfer  functions 
of  the  system  for  the  original  uncoupled  roll  channel,  Equations  (IV.B.3-3) 
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through  (IV.B.3-9),  the  following  Equations  were  derived  and  can  he 
modeled  In  a  ten  order  system  for  the  roll  channel: 


1*1  =  8.034626ft 

(IV.C4-4) 

9  =P 

(IV.C.4-5) 

5k  =  -17.69  -8X  +  1?.69c 

(IV.C.4-6) 

Y  =  -0.75325P  -0.88352<P  +0.35 165X  -5Y 
-0.403346ft  -+0.8835290 

(IV.C.4-7) 

X|  =  -6X,  +0.6829 1 6ft 

(IV.C.4-8) 

=  -0.1 027 IP  -0.120489  +0.04795X  ♦  14.3182Y 
-14.18184Xt  -156,  -0.140626ft  +0.120489c 

(IV.C.4-9) 

S2  =  +0.1 027 IP  +0.120489  -0.04795X  -  14.3182Y 
+14.18184X,  +156,  +0.1 40626r -0.120489c 

(IV.C.4-10) 

=  +0.10271P  +0.120489  -0.04795X  -  14.3182Y 
♦  14.18184X,  +156,  +0.140626ft -0.120489c 

(IV.C.4- 1 1 ) 

&4  =  -0.10271P  -0.120489  +0.04795X  ♦  14.3182Y 
-14.18184X,  -156,  -0.140626r +0.120489c 

(IV.C.4-12) 

£r  =  2698.836,  -2698.836a  -2698.836,  +2698.8364 
-188.46ft 

(IV.C.4- 1 3) 

We  thus  have  a  multiple-input-multiple-output  (MIMO)  system 
with  the  following  orden 

-  X  :  10*1,  the  time  derivative  of  X  matrix. 

-  F  :  10*10,  the  open  loop  dynamic  matrix  or  plant. 

-  X  :  10*1,  the  variables  assignments  matrix. 

-  G  :  10*5,  the  control  distribution  matrix. 

-  U  :  5*1,  the  control  input  matrix. 
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The  above  control  and  output  matrices  or  the  Uncoupled  Modified 
Unimpaired  Roll  Channel  Autopilot  are  presented  at  the  Appendix  B. 

The  Eigenvalues  found  executing  the  OPTSYS  Control  Program, 
that  Is  the  OET(sI-F)  values,  show  that  the  open  loop  system  is  stable, 
since  the  s-plane  poles  are: 

St  =  -185.057  ♦  JO 

52  =  -9.27462  ♦  j 29.2679 

53  = -9.27462  -  j 29.2679 

54  =  -2.35821  ♦  J 2.70296 

55  = -2.35821  -  J2.70296 

S*  =  -8.89209  ♦  JO 
ST  =  -5.18517  ♦  JO 
S*  =  -15.000  ♦  JO 

S9  =  -15.000  ♦  JO 

3,o=  -15.000  ♦  JO 


Applying  the  same  step  input  function  as  in  Section  IV.C.3  which 
represents  the  '1  radian  command’,  and  for  zero  rad  deflection  command 
for  the  four  control  surfaces  (i.e,  61CJ  S2o  ho  $«c  =  0)  the  time  response 

plots  of  the  following  state-variables  of  Interest  were  obtained  and 
presented  In  the  Figure  4.12  through  4.15. 

a.  Roll  Angle  Op  or  F),  Figure  4.12,  rad  vs.  time  which  has: 

(DA  time  constant  0.6  seconds 
(2)  A  5.3%  overshoot. 


(roll  angular  rate) 
(roll  angle) 


(Input  command  In  the  actuator) 
(roll  tall  Incidence) 


(3)  Zero  steady-state  error  and  so  the  guidance  will  not  be  effected. 

(4)  Delay  time  0.48  sec  and  settling  time  2.2  sec  when  the  steady 
state  roll  angle  Is  achieved. 

b.  Roll  Angular  Rate  (P  or  p),  Figure  4.13,  rad/sec  vs.  time, 
which  has: 

(DA  maximum  body  angular  rate  1.58  rad/sec  at  0.4  second. 

(2)  A  4.2*  undershoot  and  a  2.4  sec  settling  time. 

c.  Roll  Tall  Incidence  (6r  or  DR),  Figure  4.14,  degrees  vs.  time, 

which  has: 

(1)  Initially  and  very  fast,  0.1  sec,  the  maximum  required  deflection 
occurs,  that  is  1.17  degrees.  Then  changes  direction  to  the  opposite 
deflection  angle  in  0.6  sec  to  -0.32  deg.  before  stabilizing. 

(2)  A  2.6  sec  settling  time  . 

d.  .  Deflections  (51f  S2,  83,  S4  or  D's),  Figure  4.15,  rad  vs.  time, 
which  has: 

(1)  The  actuators  8t  and  64  have  the  same  sign  (positive),  as  was 
expected  according  to  the  sign  convention  of  Figure  2.2,  and  are 
first  deflected  to  a  0.02  rad  angle  and  then  to  -0.05  radians  after 
0.6  seconds.  Finally  the  steady  state  zero  angle  of  deflection  is 
achieved  at  2.6  seconds. 

(2)  The  62  and  83,  have  the  same  response  sa  above  but  with  a  negative 
sign. 

As  was  shown  the  relative  stability  for  zero-angle-of-attack  is 
satisfactory  for  the  modified  unimpaired  system. 

A  close  observation  of  Figures  4.12  throught  4.15  and  the  above 
analysis,  Indicates  that  the  corresponding  obtained  responses  are 
Identical,  with  those  of  the  original  roll  channel,  Figures  4.9,  4.10  and 


4.11.  So,  the  performance  of  the  designed  system  meets  the  pervlously 
mentioned  requirements. 

In  table  4.2  a  summary  of  the  stability  margins  is  presented, 
obtained  from  the  frequency  response  plots  calculated  up  to  100  rad/sec. 


TABLE  4.2 


PHASE  CROSSOVER  FREQUENCIES  AND  GAIN  KARGINS; 
UNCOUPLED  MODIFIED  ROLL  CHANNEL(UNIMP AIRED); 
CIRCULAR  AIRFRAME 


PHASE  CROSSOVER 
FREQUENCY(RAD/SEC) 


GAIN 

MABSINCBB) 


ROLL  ANGLE  (?) 

ROLL  ANGULAR  RATE  (P) 
ROLL  TAIL  INCIDENCE  (6p) 


10.0461 

10.0461 

10.0461 


60.2165 

60.2165 

42.1172 


We  procceed  now  to  the  analysis  of  the  Impaired  control  surfaces 
for  the  modified  system. 


D. 


CONTROL  EFFECTIVENESS  DUE  TO  THE  CONTROL  SURFACES 
IMPAIRMENT  AND  COMPARISON 


The  time  response  plots  for  the  system  are  obtained  in  this  section 
for  two  types  of  damages.  The  first  type  Is  the  locked  or  hard  over 
actuator  at  a  small  angle  and  the  second  type  Includes  the  centered  or 
inoperative  one,  two  and  three  control  surfaces. 

Note  that  the  first  phase  of  the  design  approash  for  the  CBTT 
autopilot  was  to  design  each  channel  independently  with  all  coupling 
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Figure  4.1 1  Roll  Tail  Incidence  vs.  Time;  Uncoupled  Roll  Channel; 
Continuous  Open  Loop  System;  Circular  Airframe; 
Step  Input  9C  3  1  rad;  a*  3  0- 
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Figure  4.15  Control  Surfaces  Deflections  vs.  Time;  Uncoupled  Modified 
Unimpaired  Channel;  Continuous  Open  Loop  system; 

Step  Input :  9C  *  1  and  61C>  S2»  «jc»  «4c  =  0;  ae  =  0. 


between  channels  removed.  Sufficient  high  frequency  attennuatlon  was 
added  for  actuator  and  missile  design  so  that  the  resulting  missile  body 
angular  rates  and  control  surfaces  motion  would  represent  a  practical 
missile  design.  A  relationship  was  established  among  the  relative  speeds 
of  response  of  the  uncoupled  channels  In  order  to  meet  CBTT  requirements. 
The  cross-coupling  due  to  control  surfaces  impairment  between  the 
longitudinal  and  lateral  motion  of  the  missile  Is  Ignored  when  each  control 
channel  is  treated  Independently.  Cross-coupling  of  the  controls  is  thus 
considered  In  Chapter  V. 

1.  Hard  Over  One  Control  Surface 

One  control  surface  (5t),  was  selected  to  be  deflected  or  hard 
over,  1.e.  It  remained  at  a  fixed  angle  while  the  other  three  control 
surfaces  had  a  variable  deflection  angle. 

Using  a  step  input  1  radian  for  the  commanded  roll  angle  Of  c)  (as 

In  the  original  and  unimpaired  models);  an  angle  of  0.2  rad  for  the 
damaged  surface  (6lc);  and  a  zero  angle  of  deflection,  for  the  remaining 
other  three  (62»  ha  andS^J.  The  time  response  plots  of  the  state 
variables  of  interest  were  obtained  and  analyzed  with  the  performance 
criteria  as  defined  In  section  (III.B.4). 

a.  Roll  Angle  (<p  or  F),  Figure  4.16,  (radians  vs.  time)  has: 

(1) A  time  constant  0.6  sec  and  0.3  sec  delay  time  . 

(2)  A  constant  0.4  rad  steady  state  error. 

(3)  A  5.8X  overshoot  at  1.2  sec. 


(4)  The  roll  angle  is  stabilized  at  an  Incorrect  angle  of  1.4  rad  instead 
of  the  commanded  1  radian  due  to  the  Impairment  of  5t  surface. 


b.  Roll  Angular  Rate  (  P  ),  Figure  4.1?  (radians/sec  vs.  time) 
which  has: 

(1)  Initially  and  due  to  the  fact  that  5t  is  locked  at  a  positive  angle  of 
0.2  rad,  to  achieve  the  the  roll  command,  the  angular  rate  increases 
vein  fast  the  first  tenth  of  the  second  to  a  2.1  rad/sec  rate.  Then,  a 
pronounced  oscillation  occurs  until  the  0.4  sec  when  the  rate  begins 
to  decrease  crossing  the  zero  rate  at  1.2  sec  corresponding  to  the 
max  roll  angle  (1.5  rad). 

(2)  The  initial  Increase  of  the  rate  is  due  to  the  Induced  rolling  moment 
caused  by  the  deflection  of  the  damaged  control  surface  while  the 
other  three  are  at  their  steady  position. 

(3)  A  6.6%  undershoot  at  1.5  sec  and  settling  time  2.5  sec. 

c.  Roll  Tall  Incidence  (6r  or  DR),  Figure  4.18  (deg.  vs.  time) 

which  has: 

(1)  At  time  equal  to  0.1  sec  the  actuator  Is  commanded  to  an  angle  4.75 
degrees.  This  Is  Is  due  to  the  fact  that  6f  Is  hard  over  at  the  0.2 
rad  angle  while  the  other  three  remain  have  a  zero  Initial  condition. 
Then  the  tall  Incidence  Is  stabilized  after  a  2  second  period  during 
which  a  fluctuation  occurs. 

(2)  A  zero  steady-state  error  and  2.2  sec  settling  time. 

d.  Control  Surface  Deflections  (5f,  S2,  83,64  or  D's),  Figure 
4.19,  rad  vs.  time,  which  has: 

(1)  The  actuator  S1f  is  deflected  at  the  angle  of  0.2  rad,  from  zero  sec 
up  to  3  seconds. 

(2)  The  commanded  input  vector  now  causes  the  two  surfaces  62  and  63 
to  turn  in  a  positive  angle  of  0.13  rad,  instead  of  a  negative  angle 
as  when  no  Impairment  existed. 

(3)  A  1.6  seconds  settling  time  to  achieve  this  angle. 


(4)  The  fourth  surface  S4  Is  commanded  to  turn  at  an  angle  of  equal 
amount  but  with  an  opposite  sign,  for  the  same  reasons. 

A  comparison  between  Figure  4.19  and  4.15  of  the  unimpaired 
system,  shows  that  to  achieve  the  roll  response  the  Impaired  sustem 
requires  more  control  surfaces  deflection  and  more  time. 

2.  Damage  One  Control  Surface 

The  Sf  control  surface  was  selected  to  simulate  the  inoperative 
actuator.  The  corresponding  state  variable  was  therefore  set  equal  to 
zero.  The  step  Input  of  one  radian  remained  the  same.  The  Initial  condition 
on  the  control  surface  was  set  at  zero.  This  type  of  damage  could  also 
represent  the  dedtructlon  of  the  control  surface. 

For  the  above  case  the  time  response  plots  of  the  state  variables 
of  interest  were  obtained  and  analyzed.  Later  In  Section  5,  a  comparison 
between  the  different  type  of  damages  and  the  unimpaired  system  is 
performed. 

a.  Roll  Angle  (  f  or  F  ),  Figure  4.20,  rad  vs  time,  which  has: 

(1)  A  0.6  second  time  constant. 

(2)  An  overshoot  equal  to  5.6%  and 

(3)  Zero  Steady  state  error. 

(4)  Delay  time  0.5  sec  and  2.1  sec  settling  time  at  which  the  the  roll 
response  is  achieved. 

b.  Roll  Angular  Rate  (  P  ),  Figure  4.21,  radian/sec  vs.  time, 
which  has: 


(1) A  maximum  body  angular  rate  1.58  rad/sec  at  0.4  sec. 

(2)  A  6.3%  undershoot. 


(3)  A  2.5  sec  settling  time  at  which  the  zero  steady  state  value  is 
achlved. 

c.  Roll  Tall  Incidence  (  6q  or  DR  ),  Figure  4.22,  deg.  vs.  time, 
which  has: 

(1)  At  the  very  beginning  of  the  response,  0.1  sec,  the  maximum 
deflection  occurs,  that  Is  1.1  degrees,  to  change  then  to  the 
opposite  deflection  angle  In  0.7  sec  to  -0.32  deg. 

(2)  A  zero  steady  state  error  and  settling  time  2.6  sec. 

d.  Control  surfaces  deflection  (6q  or  DR),  Figure  4.23,  deg.  vs. 

time  which  has: 

(1)  The  surface  6,  as  damaged  remains  constantly  at  zero  angle  of 
deflection. 

(2)  The  actuators  S2  and  83  are  deflected  at  a  negative  angle  0.026  rad 
and  are  stabilized  after  2.4  seconds. 

(3)  The  actuator  S4  Is  deflected  at  an  equal  but  opposite  angle  as  It  was 
expected  according  to  the  selected  sign  convention. 

A  comparison  between  the  deflection  angle  In  Figures  4.15,  4.19 
and  4.23  indicates  that  for  the  same  step  Input: 

(1)  When  one  control  surface  Is  Inoperative  or  shot  away  the  system 
commands  the  undamaged  surfaces  at  larger  angtes  (In  comparison 
with  the  unimpaired  system). 

(2)  When  one  control  surface  is  Inoperative  or  shot  away  the  system 
commands  the  undamaged  surfaces  at  smaller  angles  than  obtained 
with  the  hard  over  type  of  damage. 

(3)  The  required  time  that  the  undamaged  control  surfaces  need  to 
achieve  the  zero  steady  state  value  is  the  same  (In  comparison  with 
the  unimpaired  system). 


In  this  case  two  control  surfaces,  on  the  same  side  of  the 


missile  were  considered  damaged  (6t  and  S2).  The  two  state  variables  5, 
and  S2  were  therefore  set  equal  to  zero.  The  step  Input  of  one  radian 
remained  the  same.  The  Initial  condition  of  the  control  surfaces  was  set 
at  zero.  Then  the  time  response  plots  of  the  state  variables  were  obtained 
and  analyzed. 


a.  Roll  Angle  ($  or  F),  Figure  4.24,  rad  vs.  time,  which  has: 

(1)  Time  constant  0.59  seconds 

(2)  An  overshoot  6.8X  at  1.25  sec  and  an  undershoot  of  0.003X  at  2.2 
seconds. 

(3)  A  zero  steady-state  error  and  a  0.5  sec  delay  time. 

(4)  The  commanded  roll  angle  is  achieved  after  2.8  seconds. 

b.  Roll  Angular  Rate  (P),  Figure  4.25,  radian/sec  vs.  time, 
which  has: 

(1)  The  peak  body  angular  rate  1.61  rad/sec  is  achieved  at  0.41  sec. 

(2)  An  8.2X  undershoot  at  1.6  sec  and 

(3)  A  2.6  sec  settling  time. 

c.  Roll  Tall  Incidence  (  5c  or  DR  ),  Figure  4.26,  deg.  vs.  time, 
which  has: 

(1)  The  peak  angle  Is  1.01  degrees  at  0.1  sec  and  then  In  0.6  seconds  a 
negative  angle  (-0.35  deg)  Is  commanded. 

(2)  A  second  6.3X  overshoot  at  1.8  sec  and 

(3)  Zero  steady  state  error  (settling  time  2.6  sec). 
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d.  Control  surface  deflection  (61t  S2,  S3,  S4  or  D's),  Figure  4.27, 
rad  vs.  time,  which  has: 

(1)  The  two  Inoperative  or  shot  away  control  surfaces  5,  and  S2  remain 
always  at  zero  angle  of  deflection. 

(2)  The  two  remaining  actuators  S3  and  S4  moving  at  opposite 
directions,  as  expected,  with  a  peak  value  of  ±0.35  radians. 

(3)  The  Impairment  of  St  and  S2  has  produced  a  more  overshoot  lags 
system  as  evidenced  In  Figure  4.27. 

(4)  The  convergent  and  oscillatory  behavior  of  the  two  unimpaired 
control  surfaces  which  attempt  to  'fill  In'  for  the  damaged  surfaces 
Is  shown  In  the  same  Figure.  Also  note  that  for  this  system  with 
two  impaired  control  surfaces  more  angle  of  deflection  is  required 
(incomparison  with  the  system  having  one  actuator  shot  away). 

4.  Damage  Three  Control  Surfaces 

In  the  last  case  three  of  the  control  surfaces  8,,  $2  and  S3  were 
damaged  or  Inoperative. 

With  a  step  input  of  one  radian  for  the  commanded  roll  angle  Op c) 
as  in  the  unimpaired  and  other  type  damaged  models,  the  plots  of  the  state 
variables  of  interest  were  obtained  and  analyzed. 

Figures  4.28  through  4.31  illustrate  the  response  of  the  output 
variables.  Figure  3.30  shows  that  the  commanded  one  radian  roll  angle  can 
still  be  achieved  In  the  absence  of  three  control  surfaces  (8t,  s2,  83). 
However,  some  oscllatlon  Is  evident  In  the  response.  The  loss  of  the 
surfaces  has  produced  a  more  lightly  damped  system  as  evidenced  In 


Figures  4.28  through  4.31.  The  oscillatory  and  convergent  behavior  of  the 
roll  rate  and  the  undamaged  surface  (8^,  shown  in  Figures  4.29  and  4.31, 
Illustrate  how  the  surface  attempt  to  ‘fill  In'  for  the  absence  of  the 
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inoperative  control  surfaces.  Specifically  the  response  of  the  roll  angle 
(Figure  4.28)  has: 

(1) A  0.55  sec  time  constant. 

(2)  A  max.  10.5*  overshoot  at  1.2  sec  and  a  0.5*  undersoot  at  2.2  sec. 

(3)  The  commanded  roll  angle  Is  achieved  after  2.8  sec. 

(4)  A  0.51  sec  delay  time . 

A  comparison  with  the  response  of  the  corresponding  deflection  of 
the  unimpaired  system  (Figure  4.15)  and  the  other  Inoperative  coses.  Figures 
4.19, 4.23  and  4.31,  Indicates  that  as  the  number  of  the  Inoperative  control 
surfaces  Increases,  the  deflection  of  the  undamaged  actuators)  also 
increases  and  more  time  is  required  to  acheive  th  steady-state  value. 


In  this  section  a  performance  comparison  of  each  type  of 
impairment  with  the  unimpaired  system,  was  performed  and  an  analysis  of 
their  corresponding  time  plots  is  presented. 

For  the  above  purposes,  a  computer  program  was  written 
(Appendix  D)  to  create  a  data  file,  for  the  state  variables  of  interest,  to 
be  used  as  input  data  in  the  use  of  the  EASYPLOT  program  of  the  NPS 
computer. 

For  the  decoupled  roll  channel  the  first  plots  are  for  the  roll 
angle  and  are  presented  in  Figures  4.32  and  4.33 

In  Figure  4.32  which  compares  the  two  types  of  impairments 
with  the  original  (unimpaired)  system  it  is  quite  obvious  that  the  worst 
case  is  the  hard  over  case  because  when  the  autopilot  is  commanded  at  a 


roll  angle  (9c)  equal  to  1  radian,  the  missile  turns  at  a  roll  angle  (9)  or 
1.4  rad  Instead  of  1  rad  as  for  the  unimpaired  autopilot. 

Other  characteristic  point  is  that  for  the  hard  over  case,  the  roll 
angle  reaches  its  final  value  of  1.4  rad  (80.2°)  in  approximately  2  seconds 
with  a  5.8%  overshoot,  without  any  steady  state  error. 

Figure  4.33  compares  the  roll  angle  of  the  unimpaired  system 
with  the  three  cases  of  control  surfaces  impairments.  The  curves  include 
the  following  four  cases: 

(1)  Original  system 

(2)  Inoperative  one  control  surface 

(3)  Inoperative  two  control  surfaces 

(4)  Inoperative  three  control  surfaces 

This  Figure  shows  that  as  the  number  of  the  damaged  control 
surfaces  increases  the  overshoot  and  the  settling  time  i crease. 

Figures  4.34  and  4.35  compare  the  roll  angular  rate  of  the 
unimpaired  system  with  four  cases  of  control  surfaces  impairments.  The 
curves  include  the  following  five  cases: 

(1)  Original  system 

(2)  Hard  over  one  control  surface 

(3)  Inoperative  one  control  surface 

(4)  Inoperative  two  control  surfaces 

(5)  Inoperative  three  control  surfaces 

In  Figure  4.34  it  is  shown  that  the  hard  over  case  induce  greater 
excursion  in  the  roll  rate  with  a  peak  value  of  2.2  rad/sec  at  0.38  sec 
being  achieved.  Due  to  the  fact  that  5,  control  surface  is  deflected  at  the 


0.2  rad  the  initial  rolling  moment  In  the  first  tenth  of  the  second 
Increases  from  zero  to  2.1  rad/sec.  The  other  control  surfaces  have  zero 
as  an  initial  condition. 

Figure  4.35  compares  the  unimpaired  system  with  three  cases  of 
control  surfaces  Impairments  (Inoperative  one,  two  and  three  actuators). 
In  this  Is  shown  that  the  Inoperative  cases  are  all  quite  similar  except  for 
the  last  case  where  only  one  control  surface  Is  available.  In  this  last  case 
a  larger  delay  time,  overshoot  and  undershoot  are  evident. 

Finally  In  Figures  4.36  and  4.3?  the  roll  tall  Incidence  is 
compared  and  Include  the  following  cases: 

(1)  Original  system 

(2)  Hard  over  one  control  surface 

(3)  Inoperative  one  control  surface 

(4)  Inoperative  two  control  surfaces 

(5)  Inoperative  three  control  surfaces 

Figure  4.36  shows  that  the  required  roll  Incidence  for  the  hard 
over  case  Is  larger  In  comparison  with  the  unimpaired  system  and  with 
the  other  three  cases  of  control  surfaces  Impairment.  In  Figure  4.37  the 
maximum  peak  overshoot  decreases  with  the  number  of  control  surfaces 
that  are  impaired. 
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Figure  4.16  Roll  Angle  vs.  Time;  Uncoupled  Modified  Roll  Channel 
with  8j  Hard  Over;  Continuous  Open  Loop  System; 
Step  Input:  9C  =  1  and  Slc=0.2  rad;  520830840=0;  ae  = 
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Figure  4.20  Roll  Angle  vs.  Time;  Uncoupled  Modified  Roll  Channel  with  5t 
Centered;  Continuous  Open  Loop  System; 

Step  Input:  9c=1  rad  and  ae  =  0. 
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Figure  4.22  Roll  Tall  Incidence  vs.  Time;  Uncoupled  Modified  Roll  Channel 
with  Centered;  Continuous  Open  Loop  System; 

Step  Input:  9c=  1  rad  and  Sio^2o^3o£4c=°;  ae  =  0* 
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Figure  4.23  Control  Surfaces  deflections  vs.  Time;  Uncoupled  Modified 

Roll  Channel  with  6j  centered;  Continuous  Open  Loop  System; 
Step  Input:  rad  and  6ic=*^2o£3os4cs0»  *e  5 
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Figure  4.24  Roll  Angle  vs.  Time;  Uncoupled  flodlfted  Roll  Channel  with  6, 
and  62  Centered;  Continuous  Open  Loop  System; 

Step  Input:  *c=l  rad  and  ae =  0. 
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Figure  4.26  Roll  Tall  Incidence  vs.  Time;  Uncoupled  Modified  Roll  Channel 
with  6,  and  52  Centered;  Continuous  Open  Loop  System; 

Step  Input:9c=l  rad  and  51c,52063c»s4c=D»  ae  =  0. 
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Figure  4.28  Roll  Angle  vs.  Time;  Uncoupled  Modified  Roll  Chonnel  with  6lP 
and  63  Centered;  Continuous  Open  Loop  System; 

Step  Input:  9c= 1  rad  and  6ic>S2c>^3o^4c=0;  ae  =  0. 
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Figure  4.29  Roll  Angular  Rate  vs.  Time;  Uncoupled  Modified  Roll  Channel 
with  and  $3  Centered;  Continuous  Open  Loop  System; 
Step  Input:  9c= 1  rad  and  51c,S2o$3o$4c=0;ae  =  0. 
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Figure  4.31  Control  Surfaces  deflections  vs.  Time;  Uncoupled  Modified 
Roll  Channel :  SU62,S2  Centered;  Continuous  OpenLoop  System; 
Step  Input:  *fc=1  rad  and  ae  =  0. 
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Figure  4.32  Comparison  of  the  Roll  Angle  between  the 

Unimpaired  System  and  the  Four  Cases  of  Impairment. 
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Figure  4.33  Roll  Angle  ;  Uncoupled  Open  Loop  Channel: 

Unimpaired  and  Damaged  Control  Surfaces. 
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Figure  4.36  Comparison  of  the  Roll  Tail  Incidence  between  the 

Unimpaired  System  and  the  Four  Cases  of  Impairment. 
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Figure  4.37  Roll  Tail  Incidence;  Uncoupled  Open  Loop  Channel: 
Unimpaired  and  Damaged  Control  Surfaces. 
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A.  INTRODUCTION 

Separate  longitudinal  and  lateral  autopilots  were  designed  in  Chapters 
III  and  IV,  for  the  Modified  Linear  Open  Loop  system. 

In  this  Chapter  the  longitudinal  and  roll  autopilots  were  combined, 
resulting  in  a  new  autopilot  capable  of  commanding  pitch,  roll  or  combined 
maneuvers. 

The  time  responses  to  an  input  vector  are  Identical  to  that  obtained 
when  the  same  input  is  put  into  the  separate  pitch  and  roll  autopilots  of  the 
unimpaired  system. 

The  response  plots  of  this  system  are  then  used  as  a  basis,  to 
investigate  the  effect  of  the  damaged  control  surfaces  on  the  combined 
system.  In  this  section  the  hard  over  type  of  damage  for  the  combined 
system  is  investigated.  Corrective  actions  are  then  developed  and  presented 
and  the  robustness  of  the  complete  system  is  finally  analyzed. 

B.  DESIGN  APPROACH,  ANALVSIS  AND  PERFORMANSE  OF  THE  SYSTEM 

For  the  design  of  the  combined  system  of  the  circular  airframe 
configuration,  the  requirements  and  assumptions  mentioned  In  the  design  of 
the  separate  autopilots  are  used.  To  achieve  the  desired  maneuver  plane 
acceleration  response  for  the  unimpaired  system,  both  pitch  and  roll 
uncoupled  autopilots  were  designed  to  have  the  desired  maneuver  plane 
speed  of  response.  In  paticular,  the  goal  for  maneuver  plane  acceleration 
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response  was  a  0.5  second  time  constant  for  the  flight  condition  of  Interest 
(i.e.,  60  Kft  altitude,  Mach  3.95.  in  the  next  phase  of  the  linear  design 
technique  the  longitudinal  and  lateral  autopilot  channels  are  coupled  using 
the  aerodynamic  model  of  Chapter  III.  Both  channels  are  cross-coupled 
through  the  control  surfaces  (l.e.,  a  command  to  the  control  surface  of  one 
channel  will  affect  the  parameters  of  the  other  and  vice  versa).  The  pitch 
and  roll  state  Equations  (III.B.5-5)  through  (III.B.5-1?)  and  (IV.C.4-4) 
through  (IV.C.4-13)  are  combined  to  yield  a  twenty-three  order  system 
which  Is  presented  In  Appendix  C. 

The  Open  Loop  dynamic  matrix  Is  a  23  x  23  matrix.  This  composite 
system  Is  formed  by  two  block  diagonal  matrices  as  Is  stated  In  Equation 
(V.B-1): 

F,  (13x13)  |  0 

F  = . I .  (VB.  1-1) 

0  |  F2  (10X10) 

That  is,  the  upper  block  composed  of  the  13x13  F  matrix  from  the 
pitch  channel  and  the  lower  block  composed  of  the  10x10  F  matrix  from  the 
roll  channel. 

The  control  distribution  matrix  6,  folowlng  the  same  pottem,  is  given 
by  o  23x10  matrix,  of  the  form  and  order 

( 6,  (13x5)  |  0 

62  (10x5) 


0 


(V.B.1-2) 


Finally  the  state  vector  X  and  the  control  Input  vector  U  matrices  are 


of  the  form  and  order. 


X  =  |  23  x  1  1 


U  =  IlOxl] 


(V.B.1-3) 


(V.B.1-4) 


Because  of  this  form,  any  input  to  the  stated  combined  system,  has  a 
time  response  identical  to  that  resulting  from  the  same  input  into  the  pitch 
and  roll  chonnels  separately.  An  example  of  this  is  shown  in  Figure  5.1 
where  the  same  input  vector  used  for  the  vertical  translation  maneuver 
longitudinally  and  the  rolling  translation  maneuver  laterally  are  combined 
and  serve  as  inputs  to  the  combined  system  resulting  in  a  response 
identical  to  those  achieved  separately  in  Figures  3.9  through  3.12  and  4.12 
through  4.15. 

The  combined  autopilot  designed  in  this  thesis  and  the  simulations 
made  to  test  the  different  step  Input  commonds  ore  accomplished  using  the 
interactive,  multivariable  design  program  called  OPTSYS. 

in  addition,  in  order  to  compare  different  control  surface  impairment 
configurations  a  computer  program  was  developed  (Appendix  D).  This 
program  creates  data  files  to  be  used  as  Input  flies  at  the  EASYPLOT 
program  by  changing  the  format  of  the  OPTPLOT  data  into  a  fixed  format 
data  file.  This  Format  is  80  characters  in  length  with  space  between  the 
values. 

Three  different  command  Inputs  are  examined  for  the  Modified 
Unimpaired  system  and  are  discussed  below. 
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Test  1:  For  this  test,  both  a  pitch  normal  acceleration  and  a  roll  angle 
were  commanded  (i.e.,  Nzc  =  1  gee,  =  1  rad  and  61c,  62c»  630  64c  s  0  rad). 

In  Figure  5.1  the  time  response  plots  of  the  normal  acceleration  and  the  roll 
angle  are  presented.  As  was  expected,  both  curves  are  identical  to  the 
corresponding  plots  obtained  when  separate  autopilots  were  used. 

Figures  5.2  through  5.4  Illustrate  the  performance  of  the  other  state 
variables.  Figure  5.4  shows  the  control  surface  deflections  for  only  two 
control  surfaces  from  each  channel.  The  remaining  four  coincide  or  are  the 
mirror  image  of  the  obtained  results.  In  the  same  Figure  the  configuration  of 
the  control  surfaces  for  each  channel  Is  provided  and  the  distinction  is  also 
provided  In  the  legend. 

Test  2:  For  this  test,  only  a  pitch  command  in  acceleration  was  applied 
( I.e.,  Nzc  =  1  gee),  while  ,  9c  s  0  rad  and  5io  62c,  63c  64c  =  0  rad.  Figures 

5.5  through  5.8  show  the  performance  of  the  system.  Observe  thot  no  roll 
angle  or  any  roll  rate  is  present,  as  was  expected,  while  a  downwards 
maneuver  is  performed  according  to  the  command  input  vector  (  downward 
maneuver  is  the  positive  one).  Figure  5.8  shows  the  response  of  the  control 
surfaces  for  both  channels.  The  control  surfaces  of  the  longitudinal  channel 
are  deflected  at  a  ±0.025  radian  angle  and  of  the  lateral  channel  remained  at 

the  initial  zero  position  while  only  a  downward  maneuver  was  commanded. 
Test  3:  In  the  lost  test,  only  0  roll  commond  is  applied  ( i.e.,  f  c  =  1  rad 

),  while  Nzc  =  0  gee  and  $tc,  S2c,  63cj  5^  =  0  rod.  Figures  5.9  through  5.12 

illustrate  the  effectiveness  of  Input  vector  on  the  state  variables  of 
interest.  In  Figure  5.12  only  the  two  control  surfaces  8,  and  S2  of  the  roll 
channel  are  plotted  while  the  corresponding  surfaces  In  the  pitch  channel 


remain  at  zero.  The  other  four  coincide  with  the  obtained  results.  The 
control  surfaces  of  the  roll  channel  are  deflected  at  a  ±0.022  radian  angle 
and  of  the  longitudinal  channel  remained  at  the  Initial  zero  position  while 
only  a  rolling  maneuver  was  commanded. 

it  should  be  mentioned  here  that  according  to  the  pole-zero  plot  for  the 
designed  combined  system,  the  combined  continuous  open  loop  system  is 
also  stable,  since  the  s-plane  poles  are: 


-159.724  +  J  18.9922 

-159.724  -J  18.9922 

-8.29043  ♦  J0.0591 1 

-8.29043  -  j 0.0591 1  (pitch  angular  rate) 

-3.75933  +  J2.51441 

-3.75933  -  J2.5I441 

-0.14393  +  J0 

-0.00005  ♦  JO  (pitch  tall  incidence) 
-0.07460  ♦  j 0.66599 

-0.07460  -  j 0.66599  (pitch  normal  acceleration) 
=  -0.14300  ♦  JO 
=  -0.14300  +  J0 
=  -0.14300  ♦  JO 

=  -0.01850  +  JO  (roll  angular  rate) 

=  -9.27462  ♦  J 29.2679  (roll  angle) 

=  -9.27462  -  J29.2679 
=  -2.35821  ♦  J2.70296 
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Figure  5.1  Normal  Acceleration  and  Roll  Angle  vs.  Time;  Combined  Pitch 
and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
Sic.  *2C  fi3c  $4C  =  o  and  Nzc  =  1  gee,  <f>c  =  1  rad;  ae  =  0. 
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Figure  S.3  Pitch  and  Roll  Tall  Incidence  ve.  Time;  Combined  Pitch 

and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs 

*10  ho  ho  $4c  =  0  and  Nzc  =  1  gee,  9C  =  1  rad;  ae  =  0. 
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Figure  5.4  Pitch  and  Roll  Deflections  vs.  Time;  Combined  Pitch 

and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs 
Sto  ^20  ^Jc»  ^4c  =  0  Njg  =  1  gee,  ’Pc =  1  rad;  a^ =  0. 
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Figure  5.6  Pitch  and  Roll  Angular  Rate  vs.  Time;  Combined  Pitch 

and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
S,*  ho  ho  $4c  =  0  and  Nzc  =  1  gee,  9C  =  0  rad;  ae  =  0. 
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Figure  5.7  Pitch  and  Roll  Tall  Incidence  vs.  Time;  Combined  Pitch 

and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
S2o  s*c»  *4c  s  0  and  Nzc  =  1  gee,  9C  =  0  rad;  ae  =  0. 
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Figure  5.9  Normal  Acceleration  and  Roll  Angle  vs.  Time;  Combined  Pitch 
and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
$2»  $30  *4C  =  0  and  Nzc  =  0  gee,  f  c  =  1  rad;  ae  =  0. 
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Figure  5.10  Pitch  and  Roll  Angular  Rate  vs.  Time;  Combined  Pitch 

and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
Sid  52o  h*  «4c  =  0  and  Nzc  8  0  900.  f  c  8  1  rad;  ae  =  0. 
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Figure  5.1 1  Pitch  and  Roll  Tail  Incidence  vs.  Time;  Combined  Pitch 

and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
ho  ho  $4c=  0  and  Nzc  =  0  gee,  <pc  *  1  rad;  ae  =  0. 
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Figure  5.12  Pitch  and  Roll  Deflections  vs.  Time;  Combined  Pitch 

and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
810  $20  *3c»  S4c=  0  and  N*c  =  0  gee,  9c  =  1  rad;  as  =  0. 


S]£  -  **2.35821 

-  J 2.70296 

(V.B-22) 

S 19  =  -8.89209 

♦JO 

(V.B-23) 

S2o  =  “5.18517 

+  jO  (roll  tail  Incidence) 

(V.B-24) 

S21  =  -15.0000 

♦JO 

(V.B-25) 

S22  -  - 1 5.0000 

♦JO 

(V.B-26) 

S22  =  -15.0000 

♦JO 

(V.B-27) 

C.  THE  EFFECTS  DUE  TO  ACTUATOR  IMPAIRMENT 

In  this  section  the  effects  on  the  system  due  to  a  damaged  actuator 
was  Investigated.  One  control  surface  (6t)  was  considered  damoged.  It  was 
locked  at  a  fixed  angle  while  the  other  three  control  surfaces  had  a  variable 
deflection  angle. 

Two  different  control  inputs  were  examined  for  this  type  of 
impairment  and  are  discussed  below.  For  each  test,  corrective  actions  were 
then  developed  and  presented. 

Test  1.  Command  Input: 

(1)  5,  for  both  channels  locked  at  0.2  rad 

(2)  82,  h*  £4  for  both  channels  have  variable  deflections 

(3)  NZc  =  1  gee 

(4)  f  c  =  0  rad 

Using  the  above  control  Input  the  time  response  plots  of  the  system 
were  obtained,  analyzed  and  presented  in  Figures  5.13  through  5.16.  In 
Figure  5.13  the  pitch  norma)  acceleration  and  the  roll  angle  are  presented. 
Commanding  a  downward  maneuver  (i.e.,  N2  =1  gee  and  f  =  0)  the  missile 


body  the  first  tenth  of  the  second  moves  downward  and  then  upward  with  its 
pitch  axis  directed  at  the  velocity  vector  until  it  reaches  the  desired 
maneuver  level.  This  Is  due  to  the  Induced  pitching  moment  effect  of  the 
locked  control  surface  and  the  undamaged  control  surfaces  (82,  h*  5 4) 
which  try  to  counteract  the  impairment  of  $).  Allthough  no  roll  motion  is 
required  to  maintain  cordlnatlon  for  this  maneuver,  as  is  shown  in  Figure 
5.13  a  roll  angle  (0.4  rad)  is  present.  In  Figure  5.14  the  pitch  and  roll 
angular  rates  are  also  presented  In  which  the  1.9  rad/sec  roll  rate  is 
achieved.  In  Figure  5.16  two  control  control  surfaces  from  each  channel  are 
presented.  The  other  two  are  the  mirror  Images  of  the  obtained  results.  The 
8f  (of  pitch  and  roll  channels)  respond  properly  to  the  command  input.  The 
82  for  the  pitch  channel  Is  deflected  0.174  radians  and  for  the  roll  channel 
0.130  radians.  Both  are  stabilized  after  1.4  seconds. 

Corrective  actions:  The  roll  effect  of  the  damaged  control  surface  on 
the  performance  of  the  combined  autopilot  was  examined  here  and 
corrective  actions  to  overcome  this  problem  were  studied  and  presented. 

Figure  5.13  shows  the  roll  angle  (?)  reaching  Its  final  value  of  0.40  rad 
(23°)  In  approximately  1.4  seconds  with  a  7%  overshoot.  The  first  try  to 
overcome  the  roll  effect  was  to  command  the  opposite  surface  (83)  in  an 
equal  angle  of  deflection  0.2  radian.  The  obtained  time  response  of  the  pitch 
normal  acceleration  and  the  roll  angle  are  presented  In  Figure  5.17.  For 
comparison.  Figure  5.17  presents  three  configurations.  It  shows  that  the 
produced  roll  angle  due  to  the  8t  control  surface  Impairment  Is  eliminated. 
But,  there  is  a  penalty  in  the  pitch  normal  acceleration  which  reaches 
approximately  4  gee  and  a  oscillatory  behavior  is  evident  before  the  steady 


state  value  of  t  gee  command  being  achieved.  So,  another  test  was 
performed  by  commanding  the  63  surface  at  an  opposite  angle  of  deflection 
(-0.2  rad).  By  this  action,  shown  in  Figure  5.18,  the  normal  acceleration 
responds  very  well  but  the  roll  angle  is  Increased  at  1.18  radians  (68°)  after 
1.75  seconds,  as  might  be  expected.  The  next  test  was  performed  using  the 
second  control  surface  (62).  In  Figure  5.19  the  time  response  of  the  same 
variables  of  the  system  are  obtained  and  presented.  The  control  surface  62 
is  commanded  at  the  same  angle  of  deflection  l.e.,  0.2  radians.  As  is 
illusrtated  in  Figure  5.19  commanding  a  downward  maneuver  the  missile 
body  reaches  the  desired  maneuver  level  without  any  roll  motion  (i.e.,  Nz  =  1 

gee  and  9  =  0). 

At  the  end  of  this  Chapter  a  transition  reconfiguration  strategy  Is 
proposed  addressing  the  topics  of  failure  detection  and  Isolation,  and 
Implementation.  That  is,  one  psslble  strategy  for  transitioning  from  the 
normal  control  laws  for  a  missile  with  all  control  surfaces  operative,  to 
one  set  of  reconfigured  control  laws  for  a  missile  with  a  failed  surface. 

Test  2.  Command  Input: 

(1) 8,  for  both  channels  locked  at  0.2  rad. 

(2)  S2,  64  for  both  channels  have  variable  deflections. 

(3)  NZc  =  Ogee. 

(4)  9C  =  1  rad. 

Using  the  above  step  input  vector  the  time  response  plots  of  the  system 
were  obtained  and  presented  In  Figures  5.20  through  5.23. 

in  Figure  5.20  the  time  response  of  the  pitch  normal  acceleration  and 
the  roll  angle  of  the  combined  system  ore  Illustrated.  Due  to  the  fact  that 
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the  two  channel  autopilots  were  cross-coupled  through  the  four  control 
surfaces  (61t  S2,  63,  64),  a  damage  in  any  actuatoKs)  will  affect  the 


performance  of  both  channels.  The  locked  control  surface  S,  at  the  one  rad 


angle  affects  the  combined  autopilot.  A  constant  Incorrect  level  was 
achlved  In  the  response  of  the  roll  angle  (l.e.,  f  =  1.4  rad  Instead  of  the 
commanded  one  rad).  In  addition,  due  to  the  Induced  pitching  moment  effect 
by  the  damaged  control  surface,  a  down  and  upward  maneuver  was 
commanded  by  the  pitch  channel  during  the  first  second  of  the  response  of 


the  normal  acceleration. 


A  change  of  the  other  state  variables  was  obvious  as  Is  shown  in 
Figures  5.21  through  5.23.  Figure  5.23,  which  represents  the  time  response 
of  two  control  surfaces  from  each  channel,  t.e.,  8t  and  52,  Indicates  the 
commanded  (simulated)  failure  of  0.2  radian  at  6t  and  the  deflection  angle 


of  0.13  radians  (7.5°)  for  the  other  three  control  surfaces.  The  control 


surfaces  not  presented  in  this  plot  are  the  mirror  Images  of  the  obtained 


ones. 


t  The  roll  and  pitch  effects  caused  by  the  damaged 


control  surface  8t  on  the  combined  autopilot  performance  was  examined  and 


corrective  actions  to  overcome  this  failure  were  studied  and  presented. 


In  Figure  5.20  the  undesirable  rolling  and  pitching  effect  caused  by 
the  damaged  6,  control  surface  Is  obvious.  The  missile  rolls  about  its 
velocity  vector  with  a  1.4  radians  roll  angle  (80°)  (Instead  of  the 
commanded  1  radian).  Also  observe  the  effect  on  the  pitch  normal 
acceleration  which  should  be  constant  at  zero  since  no  down  or  upward 


maneuver  was  commanded. 


To  overcome  these  failures  several  runs  of  the  Controls  Program 
OPTSYS  were  performed  to  find  an  appropriate  corrective  action  for  the 
control  system.  The  most  Interesting  cases  are  Illustrated  in  Figures  5.24 
and  5.25. 

In  Figure  5.24  the  roll  angle  (f)  and  the  pitch  normal  acceleration  (nz) 
were  obtained  and  presented  for  three  cases  to  compare  their  responses.  A 
comparison  was  performed  between  the  original  system,  the  impaired 
system  with  the  ^  effector  locked  at  an  angle  of  0.2  radians,  and  the 
system  with  the  S2  control  surface  commanded  at  an  angle  equal  but 
opposite  (-0.2)  to  the  defective  surface.  As  is  Illustrated  in  Figure  5.24,  if 
such  an  input  was  commanded,  the  missile  body  rolls  even  more  at  a  2.2 
radians  (126°)  angle  as  might  be  expected.  Also,  the  pitch  normal 
acceleration  increases  and  achieves  the  large  value  of  -4.2  gees  before 
stabilizing  after  1.4  seconds. 

In  Figure  5.25  the  roll  angle  (f)  and  the  pitch  normal  acceleration  (nz) 
were  obtained  and  presented  but  now  the  second  control  surface  S2  was 
commanded  to  be  deflected  at  the  same  angle  as  the  S,,  (0.2  radians).  As  is 
Illustrated  In  Figure  5.25,  both  effects,  on  the  roll  angle  and  the  pitch 
normal  acceleration,  were  eliminated  and  their  time  response  plots 
indicates  that  the  missile  body  performed  according  to  the  command  input. 

In  comparing  output  time  responses,  as  was  done  by  these  two  tests,  it 
Is  apparent  that  the  Impaired  missile  was  able  to  track  the  unimpaired 
missile  as  long  as  there  was  another  control  surface  to  provide  the 
appropriate  corrective  action  and  control  power  of  the  failed  surface. 


In  the  case  where  one  or  more  control  surfaces  are  inoperative  or  shot 
away  and  are  not  considered  here,  similar  effects  are  expected.  This  is  due 
to  the  fact  that  as  far  as  all  the  control  surfaces  are  operative,  are 
commanded  by  the  system  such  that  one  will  counteract  the  induced  moment 
effect  of  the  other.  But  when  one  control  surface  Is  Inoperative  or  shot 
away  undesirable  roll,  pitch  or  combined  effects  will  be  present. 

D.  PROPOSED  RECONFIGURATION  STRATEGY 

1.  introduction 

This  section  presents  one  possible  strategy  for  transitioning  from 
the  control  laws  for  a  missile  with  all  surfaces  operative,  to  a  set  of 
reconfigured  control  laws  for  an  missile  with  a  failed  surface.  The 
discussion  begins  with  a  definition  of  reconfiguration  and  why  It  is 
important.  A  transition  stategy  Is  then  proposed,  addressing  the  topics  of 
failure  detection,  isolation,  and  Implementation. 

2.  Reconfiguration 

Reconfiguration  Is  a  technique  which  shows  promise  in  restoring 
missile  stability  and  a  degraded  level  of  controllability  In  the  event  a 
control  surface  becomes  Inoperative.  This  technique  Implies  that  the  control 
laws  governing  missile  motion  compensate  for  the  failed  surface  by  using 
the  remaining  surfaces  to  maintain  control.  The  control  laws  of  Chapter  III, 
IV  and  the  corrective  actions  of  Chapter  V  provide  the  Independent  control 
of  each  surface  necessary  to  Implement  reconfiguration  which  deals  with 


The  development  of  reconfigurable  control  laws,  such  as  those 
developed  in  this  thesis,  must  include  a  method  of  failure  identification  and 
implementation  of  a  revised  control  law.  For  a  given  failure  it  thus 
necessary  to  determine  when  the  new  lows  should  be  used  and  how  to 
switch  to  them  from  from  normal  control  laws.  A  failure  detection  and 
Isolation  scheme  must  be  developed  to  determine  when  and  which  surface 
has  failed.  The  proposed  failure  detection  and  Isolation  system  would  use  a 
set  of  redundant  sensors  and  on  estimator,  to  determine  If  the  actuators  are 
still  operative.  Should  an  actuator  fall  to  move  a  control  surface,  the 
sensors  notify  a  control  surface  monitor.  The  monitor  validates  the  fall 
signal  and  determines  which  control  surface  Is  Inoperative. 

The  control  surface  monitor  is  the  heart  of  the  reconfigurable  flight 
control  system.  When  notified  of  a  failure,  the  monitor  Informs  the  flight 
control  computers  of  the  problem  and  provides  the  necessary  Information 
concerning  which  set  of  new  control  laws  (corrective  actions)  should  be 
used  to  calculate  the  appropriate  surface  deflections.  The  off-line  flight 
control  computer  Is  programmed  with  the  new  control  laws  and  then  signals 
the  on-line  system  that  It  Is  ready  to  assume  control  of  the  missile.  Thus 
the  autopilot  Is  again  ready  to  control  a  stable  missile  with  some  degree  of 
control  capability. 
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Figure  5.13  Hard  Over  5t  Type  of  Impairment: 

Normal  Acceleration  and  Roll  Angle  vs.  Time;  Combined  Pitch 
and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
6lc=0.2rad,  S2c»  63c,  64c  =  0  and  N2C  =  1  gee,  9c  =  0  rad;  ae  =  0. 
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Figure  5.14  Hard  Over  6t  Type  of  Impairment: 

Pitch  and  Roll  Angular  Rate  vs.  Time;  Combined  Pitch 
and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
«ie=o.2rad.  5ac,  5 jo  «*c  =  o  and  n2C  =  I  gee,  fc  =  i  rad;  ae  =  0. 
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Figure  5.15  Hard  Over  5t  Type  of  Impairment: 

Pitch  and  Roll  Tail  Incidence  vs.  Time;  Combined  Pitch 
and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
61c=0.2rad,  6Zc,  SZo  5*  =  0  and  Nzc  =  1  gee,  f  c  =  0  rad;  ae  =  0. 
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Figure  5.16  Hard  Over  St  Type  of  Impairment: 

Pitch  and  Roll  Deflections  vs.  Time;  Combined  Pitch 

and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 

Slc=0.2rad,  S2C,  ho  8  0  and  Nzc  =  1  gee,  =  0  rad;  ae  =  0. 
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Figure  5.17  First  Corrective  Action  and  Comparison  of  the  Pitch 

Normal  Acceleration  (Nz)  and  the  Roll  Angle  W»  between 
Three  Different  Effectors  Configuration. 
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Figure  5.18  Second  Corrective  Action  and  Comparison  of  the  Pitch 

Normal  Acceleration  (Nz)  and  the  Roll  Angle  (f),  between 
Three  Different  Effectors  Configuration. 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


LEGEND 

NZ. ORIGINAL  SYSTEM 


a  NZ.H.OVER  D1  AND  2  =  0.2 
-  PHIt  ORIGINAL  SYSTEM 

'  ** "  PHI  tH. OVER’  D*  T  =  0.2" . 

'  O* '  PHltH.0VER  D'  V  AND'  2  =  0.2  ' 


1.2  1.6 
TIME-SEC 


Figure  5.19  Final  Corrective  Action  and  Comparison  of  the  Pitch 

Normal  Acceleration  (Nz)  and  the  Roll  Angle  (f),  between 
Three  Different  Effectors  Configuration. 
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Figure  5.20  Hard  Over  Type  of  Impairment: 

Normal  Acceleration  and  Roll  Angle  vs.  Time;  Combined  Pitch 
and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
51c=0.2rad,  S2»  $3c»  $4c  =  0  and  N2C  =  0  gee,  =  1  rad;  ae  =  0. 
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Figure  5.22  Hard  Over  8,  Type  of  Impairment: 

Pitch  and  Roll  Tall  Incidence  vs.  Time;  Combined  Pitch 
and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
_ 6ic=0.2rad,  62c>  $3c>  64c  =  0  and  Nzc  =  0  gee,  »pc  =  1  rad;  ae  =  0. 
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Figure  5.23  Hard  Over  Type  of  Impairment: 

Pitch  and  Roll  Deflections  vs.  Ttme;  Combined  Pitch 


and  Roll  Channels;  Continuous  Open  Loop  System;  Step  Inputs: 
51c=0.2rad,  52»  s4c  =  0  and  Nzc  =  0  gee,  «f»c  =  1  rad;  ae  =  0. 
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Figure  5.24  First  Corrective  Action  and  Comparison  of  the  Pitch 

Normal  Acceleration  (Nz)  and  the  Roll  Angle  (*f)»  between 
Three  Different  Effectors  Configuration. 
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Figure  5.25  Second  Corrective  Action  and  Comparison  of  the  Pitch 

Normal  Acceleration  (Nz)  and  the  Roll  Angle  (9)»  between 
Three  Different  Effectors  Configuration. 


E.  ROBUSTNESS  OF  THE  UNIMPAIRED  SYSTEM 
1.  introduction 

Up  to  this  point  the  continuous  open  loop  system  was  Investigated 
and  It  was  found  to  he  stable,  as  is  indicated  by  Equations  (V.B-5)  through 
(V.B-27).  In  this  section  the  concept  of  feedback  Is  Introduced  and  the 
robustness  of  the  system  Is  examined. 

Design  criteria  for  feedback  control  systems  can  be  lumped  into 
three  broad  categories:  performance,  stability,  and  robustness. 

Performance  criteria  are  generally  specified  as  a  requirement  for  the 
control  system  to  follow  commands  or  reject  disturbances.  The  most 
commonly  used  stability  criteria  are  stability  margins  which  generally 
refer  to  the  amount  of  gain  or  phase  variations  from  a  design  condition  a 
system  employing  feedback  control  can  experience  before  becoming 
unstable. 

In  early  feedback  control  applications,  reducing  sensitivity  was 
the  only  control  design  criteria,  particularly  reducing  sensitivity  to  noise. 
Over  the  last  several  years,  sensltlvtty  Issues  have  attracted  new  attention 
and  have  generated  a  new  requirement,  robustness.  A  control  system  Is 
robust  (Insensitive)  If  performance  Is  maintained  In  the  presence  of 
variations  in  the  plant  from  the  design  model.  These  variations  can  be 
modeling  approximations,  modeling  uncertainties,  failures  of  control 
elements  (sensors,  actuators),  noise,  non-linearities,  etc. 

These  criteria  are  not  mutually  exclusive.  Good  stability 
characteristics  In  well  designed  systems  are  usually  accompanied  by  good 
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performance.  Robustness  must  be  considered  relative  to  the  variations  from 
the  design  model  that  the  control  system  will  experience. 

2.  The  Benefits  of  Feedback  and  the  Singular  Value  Design  Technique 


d 


We  deol  with  the  stondard  feedback  configuration  Illustrated  in 
Figure  5.26.  It  consists  of  the  Interconnected  plant  (6)  and  controller  00 
forced  by  commands  (r),  measurement  noise  (n),  and  disturbances  (d). 

The  noise  In  this  case  Is  assumed  to  be  zero.  The  disturbances  are 
ossumed  to  be  reflected  to  the  measured  output  (y). 

All  signals  ore  multivariable,  in  general,  and  both  nominal 
mathematical  models  for  G  and  K  ore  finite  dimensional  linear  time 
Invariant  systems  with  transfer  functions  matrices  G(s)  and  K(s).  Then  It  Is 
well  known  that  the  configuration,  if  it  Is  stable,  has  the  following  major 
properties: 

o.  Input-output  behavior 

y  =  G  K  (I  ♦  G  K)"l  (r  -  n)  ♦  (I  ♦  G  K)'1  d  (V.E.2-1) 
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e  =  r- y 

=  (I  ♦  G  K)- 1  (r  -  d)  +  G  K  (I  ♦  6  K)“ 1  n  (V.E.2-2) 

b.  System  sensitivity 

AH  d  =  (I  ♦  G’  K)~ 1  AH  0]  (V.E.2-3) 

In  Equation  (V.E.2-3),  AH  cj  and  AH  0i  denote  changes  In  the  closed 
loop  system  and  changes  In  a  nominally  equivalent  open  loop  system, 
respectively,  caused  by  changes  In  the  plant  G,  l.e.,  G*  =  G  +  AG. 

Equations  (V.E.2-1)  through  (V.E.2-3)  summarize  the  fundamental 
benefits  and  design  objectives  Inherent  In  feedback  loops.  Specifically, 
Equation  (V.E.2-2)  shows  that  the  loop's  errors  In  the  presence  of  commands 
and  disturbances  can  be  made  "small*  by  making  the  sensitivity  operator, 
(I  ♦  G  K)"i,  "small".  Equation  (V.E.2-3)  shows  that  loop  sensitivity  Is 
improved  under  these  some  conditions,  provided  AG  Is  not  too  large 
from  G. 

With  the  Increased  Interest  in  MIMO  systems  numerous  methods  of 
design  have  been  employed  to  obtain  suitable  system  performance  and 
robustness  with  varying  degrees  of  success.  One  prlmory  method  of  design 
is  to  keep  the  plant  os  decoupled  as  possible  throughout  the  design  so  that 
each  individual  element  may  be  controlled  Independently  and  designed 
essentially  as  a  single  loop  system.  In  another  procedure  the  multiloop 
system  Is  modified  Into  a  system  that  has  diagonal  elements  that  are  much 
larger  than  any  off-diagonal  elements.  This  diagonally  dominant  system  is 

195 


then  in  a  form  where  conventional  Nyqulst  type  techniques  can  he  employed 
in  the  analysis.  A  third  common  M1MO  design  method  Is  that  of  the  Linear 
Quadratic  (LQ)  method.  This  method  uses  a  quadratic  cost  functional  and 
optimization  principles  to  allow  the  designer  to  design  for  various 
performance  levels  by  adjusting  the  matrix  weighting  terms  used  In  the 
cost  function.  The  major  difficulty  with  all  of  the  above  methods  Is  that 
they  are  not  necessarily  robust.  This  Is  especially  true  for  cross-coupling 
terms  between  loops. 

As  was  mentioned  before,  a  simple  Interpretation  of  robustness  is 
the  ability  of  the  system  to  tolerate  design  pertubatlons.  These 
perturbations  could  be  In  the  form  of  actuator  failures,  plant  parameter 
uncertainty,  unmodeled  dynamic  or  nonlinear  terms,  or  any  one  of  many 
other  perturbations  to  the  nominal  design  of  the  system. 

The  analysis  of  the  robustness  of  the  system  was  obtained  using 
the  Pole  Placement  and  Robustness  Design  Program  (POPLAR)  developed  by 
Gordon,  V.C.  [Ref.  51. 

Utilizing  the  above  program  and  using  the  data  file  shown  In 
Appendix  E  for  the  state-feedback  autopilot,  the  unimpaired  system  was 
found  to  be  robust.  The  standord  full  state-feedback  design  was  carried  out 
using  the  N.P.S.  program  OPTSYS. 

Figures  5.27  and  5.26  show  the  Minimum  Additive  Output  Singular 
Value  (SVADMO)  vs.  frequency,  and  the  Minimum  Additive  Input  Singular 
Value  (MIN.ADD.IN.SV)  vs.  frequency  respectively  for  the  state-feedback 
designed  autopilot.  It  is  noted  that  for  0.1  to  100  db  frequency  the  output 


F*6  singular  value  Is  too  low  (  0.03  rad  )  but  the  input  transfer  function 
6*F  Input  minimum  singular  value  Is  above  the  0.5  rad  region. 

Therefore  the  closed  loop  unimpaired  system  Is  robust  (at  least  in 
the  Input).  Further  analysis  would  be  needed  If  robustness  were  required  In 
both  the  Input  and  output  case  simultaneously. 

In  the  case  of  state-feedback  autopilot  and  for  a  singular  value  of 
0.5  rad  the  pole  placement  and  robustness  design  routine  places  the  poles  as 
follows: 

Ordered  (Input)  Computed  (output) 


-1706.00000 

J0.00000 

-1706.58980 

J0.00000 

-285.39990 

J0.00000 

-285.42944 

J0.00000 

-219.48999 

J0.00000 

-158.78516 

J0.00000 

-117.77000 

jo’ooooo 

-49.74738 

J0.00000 

-56.00000 

J0.00000 

-49.74738 

J0.00000 

-21.84399 

J0.00000 

-27.69325 

J0.00000 

-7.05800 

-J  1.59000 

-27.59982 

J0.00000 

-7.05800 

j  1.59000 

-8.57168 

-J2.58016 

-4.97000 

J0.00000 

-8.57168 

J2.58016 

-4.97000 

J 0.00000 

-6.42040 

J0.00000 

-4.88000 

J 2.59000 

-6.35900-J5.27725 

-4.88000 

J  2.59000 

-6.35900 

J5.27725 

-2.90000 

J0.00000 

-5.99561 

J0.00000 

-0.77000 

J0.00000 

-4.27422 

J0.00000 

-0.47000 

J0.00000 

-2.33516 

J0.00000 

-0.47000 

J0.00000 

-0.98429 

J0.00000 

-0.47000  J0.00000 

-0.40000  J0.00000 

-0.10000  jO.OOOOO 
-0.07000  -j 6.60000 
-0.07000  j 6.60000 
-0.00001  jO.OOOOO 
-4.97000  jO.OOOOO 


-0.47444  jO.OOOOO 
-0.47427  jO.OOOOO 
-0.47415  jO.OOOOO 
-0.14287  jO.OOOOO 
-0.07482  -j 6.66396 
-0.07482  j6.66396 

-0.00153  jO.OOOOO 


198 


V*  V*  > 


>  <! 


1 


kr 


vrl 


M 


vl 


LEGEND 
MIN  ADD  IN  SV 


0.0  10.0  20.0  30.0  40.0  50.0  GO.O  70.0  80.0  90.0  100.0 

FREQUENCY 


Figure  5.28  MIN  ADO  IN  SV  vs.  Frequency;  Copied  Pitch  and  Roll  Channel 
Autopilot;  State-Feedback  Design;  Continuous  Closed 
Loop  Unimpaired  System;  Ctrcular  Alrfram. 


mm 


A.  CONCLUSIONS 

The  purpose  of  this  thesis  was  to  examine  the  control  deflclecles  In  a 
Benk-to-Tum  (B.T.T.)  cruise  missile,  due  to  specific  Impairments  at  its 
control  surfaces  during  the  Intermediate  phase  of  Its  mission.  The  following 
are  the  principal  conclusions  based  on  this  worK. 

(1)  The  original  and  modified  classical  designed  autopilots  which  involve 
separate  inputs  for  each  control  surface  were  proved  to  have 
identical  performance  for  the  longitudinal  and  lateral  channels. 

(2)  The  investigation  of  specific  type  of  damages  at  the  B.T.T.  missile  s 
control  surfaces,  showed  that  the  hard  over  case  ( l.e.,  one  control 
surface  locked  at  a  certain  deflection  angle)  shows  greater  excursion 
In  the  normal  acceleration  and  roll  angle  with  large  values  being 
achieved,  than  the  other  type  of  impairments. 

(3)  It  was  showed  that  the  Impaired  missile  was  able  to  track  the 
unimpaired  missile  as  long  as  there  was  another  control  surface  to 
provide  the  appropriate  corrective  action  and  the  control  power  of  the 
failed  surface. 

(4)  A  discussion  was  given  of  a  possible  method  for  transitioning  from 
normal  control  laws  to  a  set  of  reconfigured  control  laws  for  a 
missile  with  failed  surface.  The  strategy  requires  the  system  to 
detect  and  isolate  a  surface  failure  and  then  transfer  control  to  an 
off-line  flight  control  computer  which  has  been  programmed  with  the 
new  control  laws  (corrective  actions). 

(5)  The  state-feedback  autopilot  was  Introduced  as  additional  dynamic 
design  in  order  to  implement  control  to  the  unimpaired  system.  This 
design,  analyzed  in  terms  of  its  transient  response,  was  found  to 
meet  the  desired  requirements. 


(6)  The  performance  of  the  coupled  pitch  and  roll  channel  was  found  to  be 
satisfactory  and  the  unimpaired  system  proved  to  be  robust. 


B.  RECOMMENDATIONS 

Finally  the  following  recommendations  should  be  taken  into 
consideration: 

(1)  More  complex  and  other  type  of  Impairments  could  be  investigated 
using  the  designed  autopilot,  in  order  to  examine  the  overall 
controllability  of  the  missile. 

(2)  In  order  to  make  the  system  simpler,  some  of  the  returning  gain  loops 
of  the  state-feedback  design  must  be  eliminated.  A  further 
investigation  must  then  be  conducted  in  order  to  examine  if  the 
performance  of  the  resulting  design  remains  robust. 

(3)  Further  analysis  should  be  conducted  if  robustness  were  required  in 
both  the  input  and  output  case  simultaneously. 


APPENDIX  A 
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PLANT  SYSTEM  AND  INPUT  MATRICES  OF  CONTINUOUS  OPEN  LOOP 
PITCH  CHANNEL  AUTOPILOT  (UNIMPAIRED  SYSTEM) 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


THE  OPTIMAL  FEEDBACK  GAIN  CONTROL  MATRIX. . .C=BINV*GT*S. . . 
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C-LOOP  OPTIMAL  REG.  E-VALUES ... DET( SI-F+GXC) . . 
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APPENDIX  B 


PLANT  SYSTEM  AND  INPUT  MATRICES  OF  CONTINUOUS  OPEN  LOOP 
ROLL  CHANNEL  AUTOPILOT  (UNIMPAIRED  SYSTEM) 
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OPEN  LOOP 
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THE  CONTROL  DISTRIBUTION  MATRIX  "G"-MATRIX  . . . 
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THE  OPTIMAL  FEEDBACK  GAIN  CONTROL  MATRIX. . .C=BINV*GT*S. . . 
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APPENDIX  C 


PLANT  SYSTEM  AND  INPUT  MATRICES  OF  CONTINUOUS  OPEN  LOOP 
COUPLED  PITCH  AND  ROLL  CHANNELS  AUTOPILOTS 
(UNIMPAIRED  SYSTEM) 
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0.00000*00 

0.00000*00 

0.00000*00 

9.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.23040*00 

0.11150*01 

0.17030*00 

-0.11150*01 

-0.17030*00 

0.00000*00 

0.00000*00 

0.10090*02 

0.00000*00 

0.00000*00 

0.14300*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00030*00 

0.23040*00 

-tf.il 150*01 

-0.17030*00 

0.11150*01 

0.17030*00 

0.00000*00 

0.00000*00 

-0.10090*02 

0.00000*00 

0.00000*00 

0.00000*00 

-0.14300*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

o.ooood*oo 

0.00000*00 

0.00000*00 

0.00000*00 

0 . 00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

-0.10000*01 

0.99710*03 

0.44320*03 

-0.44320*03 

0.44320*03 

-0.44320*03 

0.00000*00 

9.00000*00 

0.00000*00 

0.00000*09 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*09 

0.00000*00 

9.00000*00 

0.00000*00 

0.00000*00 

0.09000*00 

0.09000*00 

0.00000*00 

0.00000*00 

0.00900*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.90350*01 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

9.10000*01 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.09000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

-0.17400*02 

-0.90000*01 

0.00000*00 

0.00000*00 

0.00090*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00900*00 

0.00000*90 

0.00000*00 

9.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

-0.75330*00 

-0.99350*00 

0.35170*00 

-0.50000*01 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.90000*00 

-0.40330*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*09 

0.00000*09 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*09 

9.00000*00 

-0.40000*01 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.49290*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.09000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

-0.10270*00 

-0.12050*00 

0.47950-01 

0.14320*02 

-0.14190*02 

-0.15000*02 

0.00000*00 

0.00000*00 

0.00000*00 

-0.14070*00 

0.00000*00 

9.00000*00 

9.00000*00 

0.00000*00 

0.09000*00 

0.00090*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.10270*00 

0.12050*00 

-0.4795D-01 

-0.14320*02 

0.14190*02 

0.00000*00 

-0.15000*02 

0.00000*00 

0.00000*00 

0.14070*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*09 

0.00000*00 

0.00000*00 

9.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.10270*00 

0.12050*00 

-0.47950*01 

-0.14320*02 

0.14190*02 

9.00000*00 

0.00000*00 

0.15000*02 

0.00000*00 

0.14070*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

-0.15000*02 

0.00000*00 

-0.14070*00 

-0.10270*00 

-0.12050*00 

0. 47950-01 

0.14320*02 

-0.14190*02 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.20990*04 

0.00000*00 

0.00000*00 

0.24990*04 

0.00000*00 

0.00000*00 

-0.19940*03 

0.00000*00 

0.00000*00 

0.00000*09 

0.00000*00 

0.00000*00 

0.00000*00 

9.00000*00 

9.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.24990*04 

0.00000*00 

-0.24990*04 
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COMPUTER  PROGRAMS 

In  order  to  plot  the  state  variables  of  Interest  and  compare  them,  the 
following  computer  programs  were  developed  In  WATFIV  and  used.  In  order 
to  change  the  format  of  the  OPTPLOT  DATA,  into  a  fixed  format  data  file,  80 
characters  in  length  with  space  between  the  values,  to  be  used  as  input  file 
for  the  EASYPLOT  program. 


PJjTgPfji  n*  I  i 


REAL  A,B,C,D 
N=1 

WHILE  (N.LE.100)  DO 

READC5.10)A 

READ(5,15)B 

READC5, 20 )C 

READ(5.25)D 

WRITE(6,30)A,B»C,D 

FORMAT (El 4. 7) 

FORMAT ( 56X ,  El  4 . 7  ) 
FORMAT  C42X, E14 . 7 ) 
FORMAT! 42X, E14 . 7 ) 
F0RMAT(8F15 . 8 ) 

N=N+1 
END  WHILE 
STOP 
END 


and 


♦ENTRY 


REAL  A,B,C,D,E 
N  =  1 

WHILE  (N.LE.100)  DO 
READC5, 10)A 
READC5, 15)B,C 
READ! 5, 20 )D 
READ(5,25)E 
WRITE(6,30)A,B,C,E 
10  FORMAT ( E14 . 7 ) 

15  FORMAT (14X,2E14.7) 
20  FORMAT ( E14 . 7 ) 

25  FORMAT ( El  4 . 7  ) 

30  F0RMAT(4F15 . 8  ) 

N=N+1 
END  WHILE 
STOP 
END 

♦ENTRY 


[■PPT*]  FT'  I  y  IR  mil*  I  1 M  iT-TUiF]  ETTiT*  K  Pm  JK  iM-TiT'  I 


REAL  A,B»C.D,E,F,G 
N=1 

WHILE  (N.LE.IOO)  DO 

READ(5,10)A 

READC5, 15)B 

READ!5,20)C 

READ! 5, 25)D 

READ(5,30)E 

READ(5,35)F 

READ!5,40)G 

WRITE!6,45)A,E,F,G 

FORMAT ( E14 . 7 ) 

FORMAT ( E14 . 7 ) 

FORMAT ( E14 . 7  ) 
FORMAT! E14 . 7  ) 
FORMAT! 56X#  E14 .7) 
FORMAT! E14. 7) 

FORMAT ! 42X, E14 . 7 ) 
FORMAT!4F15.8) 

N=N+1 
END  WHILE 
STOP 
END 


GENTRY 


»' ?  ’77t*  vi» 


♦  JOB 

i 

REAL  A,B,C,D,E,L,G 

2 

N  =  1 

3 

WHILE  !N . LE. 100)  DO 

4 

READ!5,10)A 

5 

READ! 5, 15)B 

6 

READ! 5,20)C 

7 

READ!5,25)D 

1 

8 

READ!5,30)E 

„  1 

9 

READ! 5, 35)L 

( 

10 

READ! 5, 40)G 

| 

11 

WRITE!6,45)A,C,D.E 

1 

1 

12 

10 

FORMAT! E14. 7) 

1 

13 

15 

FORMAT! E14. 7) 

1 

J 

0 

14 

20 

FORMAT !56X«E14.7) 

! 

15 

25 

FORMAT! 42X, El 4 . 7 ) 

a 

16 

30 

FORMAT !42X,E14.7) 

f 

L 

/ 

17 

35 

FORMAT! El  4. 7) 

t 

l| 

18 

40 

FORMAT ! E14 . 7 ) 

19 

45 

FORMAT! 4F15.8) 

£ 

20 

N  =  N+1 

U 

21 

END  WHILE 

E 

22 

STOP 

i 

B 

23 

END 

) 

& 

a 

> 

♦ENTRY 

I 

t 

| 

3 

(<•«' 

!«?/ 

|'*,V 


g: 


! 

tl1! 


1$ 

v* 

y'.i 


# 


$ 

$ 


ip 

»‘rf' 

#•* 


8ft 


m, 

>*.!»< 
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PV 
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APPENDIX  E 


COUPLED  STATE-FEEDBACK  DESIGN  INPUT  DATA 
FOR  POPLAR  PROGRAM 
(UNIMPAIRED  SYSTEM) 

I  1 

100.0  0.1  0.1 
1.0  0.00  0.0 

0.5  0.0  .1  0 

0  10  0  S  1  3  1000  0 

232323102323102323 

-0. 15000000000000*03  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.15000000000000*03  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.13531000000000*01 
0.00000000000000*00  0.10000000000000*01  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  &. 00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00-0.05720000000000*01  0.00000000000000*00 
-0.50000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000*00000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000009*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.10000000000000*01  0.00000000000000*00  0.00000000000000*00 
-0.50250000000000*02  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00. 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0. 00000000000000*00-0. 04150000000000*02 
-0. 15920000000000*00  0.00000000000000*00  0.00000000000000*00  0. 17007000000000*01 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0. 10000000000000*01-0.70190000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0. 00000000000000*00-0. 4*3*0000000000*02 
-0.14920000000000*00-0.10*33*00000000*04  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000009000*00  0.00000000000000*00  0.00000000000000*00 
-0. 10940000000000*03-0. 9*200000000000*02  0.94200000000000*02-9.94200000000000*02 
0.94200000000000*02  0.00000000000000*90  0.00090000000000*00  0.00000000000000*90 
0.00000000000000*00  0.00000000000000*09  0.00000000009000*00  0.00009000000000*00 
0.00000000000900*00  0.00000000900000*00  0.00000000000000*90  0.00000000000000*00 
0.23042000000000*00  0.11 154000000000*01  0 . 1 7030000000000*00-0 . 1 1 154000000000*01 
-0.17030000000000*00  0.00900000000000*00  0.00000000900000*00  0.10092000000000*02 
-0.14300000000000*00  0.00000000000000*09  0.00000000000000 


101 


ID' 


0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00009000000000*00 
0.00000000000000*00  0.00000900000000*00  0.00000000000000*00-0.23043000000000*00 
-0.11154000000000*01-0. 17030000000000*00  0.11154000000000*01  0.17030000000000*00 
0.00000000000600*00  0.00000000000000*00-0 . 10052000000000*02  0.00900000000000*00 
-0.14300000000000*00  0.00000000000000*00  0.00900000000000*00  0.00090000000000*00 
0.00009000000000*00  0.90000000000000*00  0.00000000000090*00  0.00900000000000*00 
0.00000000009000*00  0.00000000000000*00  0.00000000000600*00  0.00000900000000*00 
9.00000000000000*00  0.00000000090000*00  0.23042000000000*00  0.11154000000000*01 
0. 17030000000000*00-0.11154000000000*01-0.1 7030000000000*00  0.90000000000000*00 
0.00000000000000*09  0.10052000000000*02  0.00000000000000*09  0.00009000000000*00 
-0.14300000000000*00  0.09000000000000*00  0.00009000000000*00  0.00009000000000*00 
0.00000000000000*09  0.00000090000000*00  0.00000900000000*00  9.00000000000000*00 
0.00000000090000*00  0.00000009000000*00  0.00000000000000*90  0.00000009000000*00 
0 . 00000000000090*00-0 . 23042000000000*00-0 . 1 1 154000000000*0 1-0.1 7010900000000*00 
0.11154000000000*01  0.17030000000000*00  0.00000090000000*00  0.00900000000000*00 
-0.10052000000000*02  0.00009000000000*00  0.00009000000000*00  9.00000000000000*00 
-0.14300000000000*00  0.00000000900000*00  0.00090000000000*00  0.00090000000000*00 
0.00000000000000*00  0.00009000000000*00  0.00009000000000*00  0.00000009000000*00 
0.00000000009000*00  0.00000000009000*00  9.00000000000000*00  0.00000009000000*00 
0.00009000000000*00  6.00000000000000*00  0.00000000000000*09  0.00009000000000*00 
0.00000090000000*00  0.00000000000000*00-0.10000000000000*01  0.00712400000000*03 
0.44321100000000*01-0.44321100000000*03  0.44321100000000*01-0.44321100000000*03 
0.00000900000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000900000000*00  0.00000090000000*00  0.00000000000000*00  0.90000000000000*00 
0.00000000900000*00  0.00009000000000*00  0.00000000000000*09  0.00090000000000*00 
0.90000000000000*00  0.00000900000000*00  0.90000000000000*00  0.00000000000000*09 
0.00000000000000*00  0.00000000000000*90  0.00000000000000*09  0.00000090000000*00 
0.00000000000000*09  0.00000000000000*09  0.09000000000000*00  0.00000000000000*00 
0.90000000000000*00  0.00600000000000*00  0.00090000000000*00  9.00000000000000*00 
9.00000000000000*00  0.000*0000000000*00  0.00000000009000*00  0.00000090000000*00 
0.00000000000000*90  0.00344200000000*01  0.00090000000000*00  0. 00000000000000*00 
0.00000000000007*00  9.00000000000000*00  0.00000000000000*09  0.00000000000000*00 
0.00000000000000*90  0.00000000900000*00  0.00000900000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000900000000*00  0.00000900000000*00  0.19000000000000*01 
0.00000000090000*00  0.00000000000900*00  0.00000900000000*00  0.00000000900000*00 
0.00000000000000*00  0.00009000000000*00  0.00000000000000*90  0.00090000000000*00 
9.00000000000000*00  0.00000000000000*06  0.00000000000000*09  0.00009000000000*00 
0.00000000000000*00  0.00000000000000*09  0.00000090000000*00  0.00000000090000*00 
9.00000000000000*00  0.00000000900000*00  0.00000000009000*00  0.00000000000000*00 
0 . 00000000000000*00  0 . 90000000000000*00  0 . 00000000000900*00-0 . 1 7400000000000*02 
-0.00000000000000*01  9.00000000000000*00  0.00000000000000*09  0.00900000000000*00 
0.00000000000090*00  0.00000000900000*00  0.00900000000000*00  0.00000000000000*00 
0.00000000009000*00  0.90000000000000*00  0.00000000900000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  9.00000000000000*00  9.00000000000000*00 
0.00000000000000*00  0.00000090000000*00  0.00000000000000*90  0.00000000000000*00 
0.00000000000000*00-0. 75325000000000*00-0. 90352000000000*00  0.35145000000000*00 
-0.50000000000000*01  0.00000090000000*00  0.00000000000000*00  0.00000090000000*00 
9.00000000000000*00  0.00090000000000*00-0.40334000090000*00  0.09000000000000*00 
0.00900000000000*00  9.00000000000000*00  0.00000000000090*00  0.00000900000000*00 
0.00000000000000*00  0.00000000000000*00  0.00090000000000*00  0.00900000000000*00 
0.00000900000000*00  0.00000000000000*09  0.00000000900000*00  0.00000009000000*00 
0.00000000000900*00  0.00000090000000*00  9.00000000000000*00  0.00000000000000*90 
-0.40000000000000*01  0.00090000000000*00  0.00000090000000*00  0.00000000000000*09 
0.00000000000000*00  0.49251000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00090000000000*00  0.00000000000000*00  0.00090000000000*00 
0.00000000000000*00  0.00000000000000*09  0.00000900000000*00  0.00009000000000*00 
0 . 00000000000000*00  0 . 90000000000000*00  0 . 00000000000000*00-9 .10271 000000000*00 
-0. 12049000000000*00  0.47550000000000-01  0. 14319290000000*02-0. 14191940000000*02 
-0.15000000000000*02  0.00000009000000*00  0.00000000000000*09  0.09000000000000*00 
-0.14047000000000*00  0.90000000000000*00  9.00000000000000*00  0.09000000000000*00 
0.00000000000090*00  0.00000000000900*00  0.00000090000000*00  0.00000009000000*00 
0.00000000000900*00  0,00000000000000*09  0.00000900000000*00  9.00009000000000*00 
0.09000000000000*00  9.00000000000000*00  0.10271000000900*00  0.12949000000000*00 
-0.47550000000000-01-0.14319200000000*02  9.14191940000000*02  0.90000000000000*00 
-0.15000000000000*02  0.90000000000000*00  0.00000000000090*00  0.1 404 7000000000* 00 
0.00000000000000*00  6.00900000000000*00  0.09000000000000*00  0.00000000900000*00 
0.00000000000000*09  9.00000000000000*00  0.00000000000000*09  0.00000000000000*00 
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o . oooooooooooood*oo  0 
0.00000000000000*00  0 
-0. 14311200000000*02  0 
-0.15000000000000*02  0 
0.00000000000000*00  0 
0.00000000000000*00  0 
0.00000000000000*00  0 
-0. 10271000000000*00-0 
-0. 14101840000000*02  0 
-0. 15000000000000*02-0 
0.00000000000000*00  0 


00000000000000*00  0. 
10271000000000*00  0 
16 10 J04 00 00000*02  0 
00000000000000*00  0. 
00000000000000*00  0. 
00000000000000*00  0. 
00000000000000*00  0. 
12069000000000*00  0. 
00000000000000*00  0. 
14067000000000*00  0. 
00000000000000*00  0. 


0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.26909300000000*04-0. 
-0.11040000000000*03 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  o. 

0.00000000000000*09  0. 
0.00000000000000*00  0. 
0.00090000000000*00  0, 
0.65720000000000*01  0. 
0.00000000000000*00  0. 
0.00000000000000*90  0. 
0.00000000000900*00  0. 
0.00000000000090*00  0. 
0.00000009000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
-0 . 94200000000000*02  0. 
0.00000000000000*00  0, 
-0.14399099909990*99  9. 
-0.23042000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
-0.23042000000000*00  0. 
9.00000000000000*00  0. 
0 . 00000000000000*00-0 . 
0.00000000000000*00  0. 
0.443211 00000000*03-0 . 
0.00000000000000*00  0. 
0.00000000000090*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00900000000000*00  0. 
0.09000000000000*00  0. 
0.00000000000000*00  0. 
0.90000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000090000*00  0.1 
0-00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00  0. 
0.00000000000000*00-0. 
0.00000000000900*00  0. 
0.00000000900000*00  0. 
-0.15000000000000*0 2  0. 
9.00000000000000*00  o.i 


00000000000000*09  0 
00000000000000*00  0 
26911300000000*04-0 


00000000000000*00  0. 
12048000000000*00-0. 
00000000000000*00  o. 
16067000000000*00  0. 
00000000000000*00  0. 
00000000000000*00  0. 
00000000000000*00  0. 
47650000000000-01  0. 
00000000000000*00  0. 
00009000000000*00  0. 
00000000000000*00  0. 
00000000000000*00  0. 
00000000009000*00  0. 
00000000000000*00  0. 
26 Ml 30000 00 00* 04  0. 
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00000000000000*90  0 
00000000000000*00  0 
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00000009000000*00  0 
00000000000000*00-0 
64200000000000*02  0 


. 00000000000900*00  0.00000000000000*00 
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.00000000000000*00  0 
.00000000000000*00  0 
.00000000000000*00  0 
.00000000000000*00  0 
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00000000000000*00  0 
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11352000000000*00  0 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


0.00000000000000*08  0.00000000000000*00  0.00000000000000*00-0.15000000000000*02 
0. 00000000000000*00-0. 120*0000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000800*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00-0.15000000000000*02  0. 12068000000000*00 
0.00000000800000*00  0.00000000000000*00  0.00800000000000*00  0.00000000000000*00 
0.00000090000000*00  0. 20*00500000000*04-0. 2t9S930000000D*0«-0.2499t30000000D*0« 
0 . 2 i 9003 00000000*04  0.80000000000000*00 

0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
O.OOOOOOOOOOOOOD*OQ  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00080000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*80  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00009000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.09000000000000*00 
0.08000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000800000000*00  9.00000000000000*00  0.00000000000000*08  0.00000000000000*00 
0.10000000000000*01  9.00000000000000*90  0.00000000000000*00  0.00000000000000*00 
0.90000000000900*00  9.00000000000000*00  0.00000000000080*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000800000000*00  0.00000000000000*90  0.00000000000000*00  0.00000000000000*00 
9.00000000000000*00  0.00000000000000*00  0.00000000000080*00  0.00000000000000*00 
0.00000000000000*08  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000080*01  0.00000000090000*00  0.00000000000000*00  0.00000009000000*00 
0.00000000000000*00  0.90000000000000*00  0.00000000000000*80  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00080000000000*00  0.00000000000000*00 
0.00000000800000*00  0.00000000000000*00  9.00000000000000*00  0.00090000000000*00 
0.00000000000000*00  0.00000000000000*90  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*80  0.00000000000000*80 
9.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.08000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.90000000000000*00  0.00000090000000*00  0.00000000000000*00 
9.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000008000000*00  0.90000000000000*00  0.00000000000000*00  0.00000000900000*00 
0.00008000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
8.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000800000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00800000000000*00  0.00000000000900*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  9.00000000000000*00  9.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  9.00000000000000*00  0.00000000000000*09  0.00000000000000*00 
0.00000000000000*00  0. OOOOOOOOOOOOOD*09  9.00000000000000*00  0.00009000000000*00 
0.10000800000000*01  0.00000900000000*00  0.00000000000000*00  0.00000900000000*00 
0.00800000000000*00  0.00000000000000*00  0.00000000000000*09  0.00000000000000*00 
0.80000000000000*00  0.00900009000000*00  0.00000000090000*00  0.09009009000000*00 
0.00000000000000*00  0.00000000000080*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  9.00000000000000*00  9.00000000000000*00  0.00000000000000*90 
0.00800000000000*00  0.00000000000000*00  0,00000000000000*09  0.00000000000090*00 
9. 19909099009000*91  9.00009009000000*90  9,00000000000000*00  0.00090000000900*00 
0.00000800000000*00  0.00000000000000*00  0,00000000000000*00  0.00009000000000*90 
0.00000000000000*09  0.00000000000000*00  0.00000000000000*00  0.00000000000000*09 
0.00000000000000*09  0.00000000000000*00  0,00000000000000*00  0.00000000000000*00 
w. 00000000000000*00  0.00000000000090*00  0,00000000000000*00  0.00000000000000*00 
0.00000900000000*00  0.00000000000000*00  0.00000090000000*00  0.00000000000000*00 
0.10000000000000*01  0.00009000000000*00  0.00000000000000*00  0.00080000000000*00 
0.00000000000000*00  0.00000000000090*00  0.00000900000000*00  0.90000000000000*00 
0.08000000000000*00  0.00000000000000*90  0.00000000000000*00  0.00000008000000*00 
0.00009000000000*00  0.00000000000000*00  0.09000000000000*08  0.90000000000000*00 
0.00000000000000*00  8.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000900000*00  9.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
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REPRODUCED  AT  GOVERNMENT  EXPENSE 


0.10000000000000*01  0 . 00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0 . 1 0000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  o.oooooooooooood*oo  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  o.oooooqoooooood*oo  0.00000000000000*00 
0. 10000000000000*01a  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.06000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.0000000000000D*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.0000000000000D*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  o.oooooooooooood*oo  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.60000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*60  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  o.oooooooooooood*oo 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.0000000000000 0*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  o.oooooooooooood*oo 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*60  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 


1*0 


CED  AT  GOVERNMENT  EXPENSE 


0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  O.OOOOOO00O0OOOD*QO  0.00000000000000*00 
0.10000000000000*01 

0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0 .00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 


0 . 00000000000000*00  0. 00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 


0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*09  9.00000000000000*00 
0.00900000000000*00  0.00000000000000*90  9.00000000000000*00  0.00900000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*09  9.00000000000000*00 
0.90000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00009000000000*00 
0.00000090000000*00  0.00000000000000*09  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000000*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  0.00000000000090*00 
0.00000000000000*00  0.00000000000000*00  0.00000000000000*00  9.00000000000000*00 
0.10000000000000*01  0.00000000000000*00  0.00000000000000*00  0.00900000000000*00 
0 . 00000000000000*00  0.00000000000000*00 
10000000000000000000000 
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